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Final report for NASA grant NAG5-1554 

This serves as a final report for grant NAG5-1554 entitled Energization and Transport of 
Ions of Ionospheric Origin in the Terrestrial Magnetosphere. The work has been predominantly 
focused on ion outflows from two data sets: 1) Prognoz 7 and Dynamics Explorer. 

The Prognoz analysis has been carried out in conjunction with Dr. Oleg Vaisberg of IKI. 
The study analyzed ion densities, temperatures, and flow velocities in the magnetotail. The work 
performed under this contract consisted of developing a program to load the raw data, compute 
the background subtraction of a strong sun pulse, and using the net counts to calculate the low 
order moments of the distribution function. The study confirms the results of ISEE with regard 
to the supply of plasma from the cusp as a major source of plasmasheet plasma and goes beyond 
this to discuss the use of ion velocities as a way to examine the motions of the magnetotail. The 
work to be reported in Advances in Space research is included in preprint form as Appendix A. 

The work on the DE/Retarding Ion Mass Spectrometer data set consumed the bulk of the 
effort. The work provided under this grant is separated into two categories: 1) classification of 
low-energy ion flows from high-latitudes, and 2) studies of the polar wind. The work on 
classification of ion flows has consisted mainly of support work on studies that begin before my 
departure from NASA Marshall Space Flight Center. The major publications that have resulted 
from this work are: 1) the determination of plasma fluxes from the cleft ion fountain Pollack et 
al.[ 1990] and the recent thorough analysis and characterization of DE/RIMS high-latitude ion 
outflows by Giles et al. [1995]. These reprints are included as Appendix B of this report. 

Although, the polar wind is indeed a component of the low energy ion outflow at high 
latitude, the history of extensive theoretical work over the last two decades of this thermal escape 
of hydrogen and helium as a result ot ambipolar diffusion through the heavier, oxygen-dominated 
topside ionosphere puts it in category all to itself. The analysis of the polar wind by the RIMS 
group in the late 1980’s culminated in the work of Chandler et al. [1991] (see Appendix C). 
This work reports the flux variability as a function of season, magnetic activity, etc. and was 
partially supported by this grant. In addition much effort has been expended under this grant to 
complete a follow on study of the thermal structure of the polar wind. Extensive display tools 
and analysis software have been developed and used in an attempt to carry out this thermal 
analysis. However the difficulties of performing multi-parameter, non-linear least squares fits 
has so far proved to be insurmountable. Therefore, the present work (see notes in Appendix C 
on work still in progress) uses a constrained fit scheme that combines the ion densities and flow 
velocities derived from Chandler et al. [1991] and a spacecraft potential derived from an 
empirical relation to the total ion density (utilizing the complete RIMS data set) to determine the 
remaining fit parameter, the ion temperature, via a least squares fit to the RIMS data. Reports 
of this work are anticipated in the near future. 
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ABSTRACT 

The plasma spectrometer SKS on the Prognoz 8 satellite was able to measure the very dilute 
plasma in the tail lobe at distances 15-20 Earth radii about 50% of the time. The energy of the 
ions, being in the range of several hundred electron volts, decreases as the satellite moves from 
the high-latitude magnetopause to the plasma sheet. The velocity spread of the ions is narrow 
in the direction of motion, corresponding to a temperature of several eV. The temperature in the 
transverse direction is 2 to 3 times higher. Comparison with magnetospheric convection models 
suggests that these ions originate at low altitudes in the region near the cusp and are separated 
in the magnetospheric velocity spectrometer. Observed velocity variations along the spacecraft 
trajectory suggest the existence of a characteristic velocity profile and a characteristic 
temporal/spatial scale. Comparison of simultaneous E/q measurements with mass-selective 
detector measurements ot protons and the alpha particles suggests that lobe plasma at these 
distances consists primarily of oxygen ions. The average ion flux density is about 3xl() 5 cm' 2 
s' 1 , giving estimated total flux over the lobe about 10 2 .s' 1 , comparable to cusp/cleft source. 



Appendix B 

Pollack et al reprint 
Giles et al. reprint 



JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 95, NO. All, PAGES 18.969^18,980, NOVEMBER 1, 1990 


A Survey of Upwelling Ion Event Characteristics 

C. J. Pollock. 1 M. O. Chandler, 1 T. E. Moore, 1 J. H. Waite, Jr, 2 
C. R. Chappell. 1 and D. A. Gurnett 3 

Ionospheric ion upwelling in the vicinity ot the dayside cleft has been studied, based primarily on data from 
the Dynamics Explorer 1 spacecraft. Using retarding ion mass spectrometer low-energy ion data and plasma 
wave instrument dc electric field data, bulk ion plasma parameters, including ion species density and field- 
aligned bulk velocity and flux, have been derived at points within a number of observed upwelling ion events for 
the ion species H . He . O . and O * + The ion species bulk parameters near the source latitude are examined 
and compared. We find that the upwelling plasma is rich in O . which typically comprises —909^ of the particle 
density, followed by H~ at somewhat less than 107r, and then He* and 6**. each comprising ~17r of 
the upwelling ion panicle density. The upwelling O* ion flux is also commonly dominant over that associated 
with the other species, with normalized values near the source region which are typically near I0 4 cm : s l . The 
fractional upward H flux is not as small as the fractional H* density due to the much larger H* upward 
flow velocities. Integration ot the product ot the normalized upward ion species flux and the upwelling ion 
occurrence probability (Lockwood et a!., 1985) over the source area yields an estimate of the source strength of 
this low-altitude cleft region magnetospheric plasma source of 2.6 x I0 :5 ions s'. 


1 . Introduction 

During recent years, it has been consistently demonstrated that 
the Earth's ionosphere acts as a significant source of magneto- 
spheric plasma [Shelley et al., 1972; Young et at .. 1977; 
Ghiehnetti et a /., 1978; Johnson . 1979; Balsiger et at .. 1980; 
GornexetaL, 1981; Lundin et at., 1982 : Collin etaL 1984; Yau 
etaL, 1984. 1985; Lennartsson and Shelley. 1986). In fact, it has 
recently been argued [Chappell et at .. 1987] that ionospheric 
plasma may fully populate the magnetosphere, with no significant 
solar wind source being required to account for observed mag- 
netospheric charged panicle populations. Recognition of the 
importance of the ionosphere as a source of magnetospheric 
plasma began with observations ot heavy ions, previously 
believed to be gravitationally bound to low geocentric altitudes, at 
magnetospheric altitudes. As such observations continued to be 
reported, in diverse magnetospheric regions such as the ring 
current [Johnson et al., 1977). the plasma sheet [Peterson et al ., 
1981 |. and the plasma sheet boundary layer and tail lobes 
| Eastman et al., 1984), it became apparent that escaping plasma 
transport from very low altitudes must be quite common. At the 
beginning of the 1980s, however, the questions of transport 
mechanisms and specific source regions were only beginning to 
be obscrvationally addressed. 

Observations of outflowing ionospheric ions have been 
reported by a number of authors in the last decade. Outflowing ion 
distributions which have apparently evolved, beginning with 
transverse (to B) acceleration and followed by magnetic mirror 
folding, have been observed near 1 R F { Sharp et al., 1977] in the 
topside ionosphere | Klumpar, 1979] and at sounding rocket alti- 
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tudes [Whalen etaL, \91&, Yau et al.. 1983; Moore et al., 1986 a; 
Kintner et al., 1989]. Since the launch of Dynamics Explorer 1 
(DE I), observation of low-energy ionospheric plasma outflows 
has become very common. The realization that the dayside mag- 
netospheric cleft region represents a significant source of outflow- 
ing plasma began early in the DE mission, based on backward 
trajectory analysis of upflowing beams observed in high-altitude 
(2-4 R f ) polar cap field lines [Waite et al., 1985). This realization 
reached maturity with the work of Lockwood et al. [1985] who 
used DE 1 retarding ion mass spectrometer (RIMS) data to study 
the statistical occurrence of various types of ionospheric outflows. 
That work identified a persistent (occurrence probability exceeds 
0.5 in source region) source of plasma emanating from the day- 
side auroral zone, showed that this source apparently favors the 
prenoon region of magnetic local time, and showed that there is a 
strong relationship between the low-latitude edge of the source 
region and the location of (usually downward) field-aligned 
currents. They labeled the observation of the distinctive RIMS 
spin-time signature associated with this outflow as an upwelling 
ion event (UWI). Moore et al. [19866] and Waite et al. [1986] 
presented case studies of an UWI, observed near 1930 UT on day 
71 of 1982. They derived bulk flow parameters for several species 
and demonstrated that the upwelling ion plasma core is hot and 
conically shaped, having undergone transverse energization at 
some lower altitude. They also noted a strong shear in the 
transverse electric field associated with the event. 

The purpose of the present study is to more broadly investigate 
the bulk properties of UWIs. using RIMS and plasma wave instru- 
ment (PWI) quasi-static electric field data. The results of the study 
yield quantification of the densities and the upward field-aligned 
velocities and fluxes, for the ion species H + , He^, 0 + , and 
+ , associated with these flows. These parameters, as measured 
near the lower-latitude edge of the events, where the flows are 
most intense, will be presented. All of the available altitude 
normalized ion flux measurements from within the upwelling ion 
events are combined with the statistical occurrence probabilities 
derived by Lockwood et al. [1985] to produce spatial source 
strength distributions. Integration over these distributions yields 
an improved measure of the total low-energy ion outflow rate 
from this source region. 
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2. INSTRUMENTATION 
2.1. The Spacecraft 

The observations presented in this report were gathered using 
instruments flown aboard the Dynamics Explorer I <DE 1 ) space- 
craft. DE I was launched on August 3. 198 1 . into a highly ellipti- 
cal polar orbit with apogee and perigee altitudes at 23,300 km 
(4.66 R h; . geocentric) and 570 km. respectively. The orientation 
of the orbit plane is fixed in the geocentric equatorial inertial 
(GEI) frame of reference, so that in the geocentric solar ecliptic 
(GSE) frame of reference, it is seen to rotate westward at a rate of 
0.9 8°/ day. This provides for approximately one full 360° rotation 
of the orbit plane in geocentric local time per year. Further, the 
DE 1 line of orbital apsides rotates about the negative orbit plane 
normal at a rate ot 120°/vear, so that the geographic latitudes of 
apogee and perigee drift retrograde at that rate. The spacecraft 
orbit, then, very nearly repeats itself, in the Earth s frame of refer- 
ence. once every 3 years. These aspects of the spacecraft orbit 
influence our ability to study UWIs, not only because the source 
region has a limited geocentric solid angular extent, but also 
because the flowing plasma seems to evolve substantially with 
altitude. Therefore, the plasma characteristics which give rise to 
the distinctive RIMS UWI signature appear to exist only in a small 
region of space sampled by the spacecraft for several months once 
every 3 years. Furthermore, the orbit forces a phase-locked 
altitude-local time sampling which may bias the occurrence statis- 
tics of the local time extent of the upwclling ion source region. 
Other aspects of the DE mission have been described elsewhere 
\ Hoffmann and Schmerling, 1981). 

2.2 Retarding Ion Mass Spectrometer 

Positive ion parameters presented in this report were obtained 
using the DE 1 retarding ion mass spectrometer (RIMS) [Chappell 
et al.. 1981]. The RIMS instrument consists of three similar 
sensor heads, one each viewing parallel and antiparallcl to the 
spacecraft spin axis (negative and positive Z heads, respectively) 
and the third viewing in the spin plane (radial head). The two Z 
heads feature circular fields of view with half angles of 45°. while 
the radial head views a field which has a 45 !J half anele in the 
plane perpendicular to the spin plane and a 20° half angle in the 
spin plane. Since the DE l spin axis is nearly perpendicular to 
the geomagnetic field, the fields of view of the two Z heads are 
centered nearly perpendicular to that field, and the field of view of 
the radial head sweeps through a nearly full range of pitch angles 
w ith the spacecraft spin. All three sensor heads were designed to 
discriminate among incoming positive ions on the basis of energy 
per charge (0—50 eV/</), using a retarding potential analyzer 
(RPA). and on the basis of mass per charge ( 1-32 aniu/</). using a 
magnetic mass spectrometer Further, each sensor utilizes two 
detectors that are sensitive to ions with mass ratios of 4 to 1 such 
that, for example, while one detector is viewing O 4 , the other is 
viewing He" . The RPA on the radial head failed in late 1981 , with 
the retarding grid having become electrically grounded to the 
spacecraft chassis so that, subsequent to the failure date, all ions 
with energies per charge greater than any positive spacecraft 
potential were passed into the magnetic mass spectrometer. The 
mass spectrometer, in turn, is characterized by a finite energy 
passband. whose width is dependent upon the species mass per 
charge and varies from approximately 250 cV for H ’ down to 16 
cV tor CT. The upwclling ion densities, field-aligned bulk 
velocities, and fluxes presented below were derived using ion data 
from the RIMS radial head. 


2.3. Quasi-Static Electric Fields 

Quasi-static electric field data presented in this report and used 
in deriving ion bulk parameters were obtained using the Z axis 
antenna ot the plasma wave and quasi-static electric field instru- 
ments [Shawhan et al . , 198 I |. I his tubular. 9-m (tip to tip) 
antenna is sensitive to the electric field component perpendicular 
to the spacecraft spin plane. The geomagnetic field is nearly 
(within 5 in 9Vf ot our cases and always within 12“ for our 
data set) contained within the DE I spin plane, so that the / 
component electric field is associated with convective plasma 
drifts which are approximately w ithin the spin plane and perpen- 
dicular to the ambient geomagnetic field. This electric field 
component is provided at a rate of 16 samples s 1 over a nominal 
range ot 0.5 m\ m 1 to 2 V rrf 1 . Characteristics of this instrument 
which are pertinent to the analysis techniques employed in this 
study will be described below. 

3. Data Analysis 

3.1. Data Selection 

The observations presented here arc taken from a subset of up- 
wclling ion events identified by M. Lockwood through a survey 
of the RIMS data as displayed in the DE I microfiche (with format 
similar to Figure 2 a displaying data for the species H * and He * ) in 
conjunction with his study of the occurrence of various tvpes of 
terrestrial ion outflows [Lockwood et al., 1985). Lockwood et al. 
used a qualitative visual signature identification technique to 
identify these upwclling ion events. The RIMS spin-time up- 
welling ion event signature has been described elsewhere (see. for 
example, Plate 3 of Lockwood et al. (I985|). This signature is 
again illustrated in the O ‘ data presented in Figure 2a. where the 
count rate in the CC channel is shown plotted in grey scale versus 
UT (along with certain orbital parameters) along the abscissa and 
the RIMS radial head viewing direction, with respect to the space- 
craft ram direction, along the ordinate. The solid and dashed lines 
running through the spectrogram from left to right indicate view- 
ing angles along the positive and negative geomagnetic field 
directions, respectively. For northern hemisphere passes, ions 
detected while viewing along the positive magnetic field direction 
are moving in the anti field direction and, therefore, away from the 
Earth. The signature ot the event is evident in Figure 2. bemnninu 
near 2320 UT and continuing until just past 2327 UT. The ramp- 
like upward flux intensification, followed by a sharp cutoff and 
return to ram-dominated flow into the instrument at the low- 
latitude edge ot the event, typifies the signature of these events 
used by Lockwood et al. in their fiehe survey. 

Having surveyed the available data obtained between October 
1981 and October 1983, Lockwood et al. identified 86 upwellinii 
ion events. These events were identifiable in the fiehe for periods 
of time ranging from 1 to 15 min, providing a total of nearly 450 
min of observation within the events. Of these 86 events, we have 
selected 39 for inclusion in the analysis for this report. This selec- 
tion was based, in part, on the requirement that the RIMS instru- 
ment be operating in a mode without aperture bias applied to the 
instrument entrance and on the further requirement that the meas- 
ured Z component of the electric field was deemed reliable. The 
RIMS instrument features an aperture bias capability, whereby a 
negative electrostatic potential may be applied to a grid at the 
entrance aperture, in order to overcome the effects of a positive 
spacecraft potential, allowing access of very low-energv ions to 
the instrument. Of the 86 upwclling ion events identified in the 
Lockwood et al. study. 20 were from times when RIMS was 
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operating in the aperture bias mode, and analysis of these events 
has been deferred pending the development of routine analysis 
tools for this mode. For another 8 events, either the digital RIMS 
data or the attitude and orbit data were unavailable. Of the remain- 
ing 58 events. 16 were rejected because the Z antenna de electric 
field data were deemed unreliable. This antenna has a length of 
9 m tip to tip, as compared with the spin plane antenna which is 
200 m long. It is known to be subject to stray fields that originate 
from differential spacecraft charging and can have large amplitude 
relative to the fields of interest. This left 42 events, from which 2 
were rejected due to a combination of no electric field data avail- 
able and highly variable, difficult to interpret ion data and I was 
rejected due to the apparent presence of a telemetry dropout dur- 
ing the event. Thus, 39 of the original 86 events, or 45%, have 
been processed to obtain ion bulk parameters along the spacecraft 
track. Of these, 5 have no electric field data, and values of E z = 0 
have been assumed. These 39 events, consisting of a total of 
nearly 200 min of observation, constitute the data set used in the 
present analysis. 


3.2. Derivation of Ion Bulk Parameters 

Due to the failure of the RIMS radial head RPA, ion energy 
spectra arc not available in the spin plane data. This situation has 
led to the development of an alternative technique for deriving ion 
bulk parameters from the RIMS radial heat data. This geometrical 
technique, developed by Chandler and Chappell 1 1 986 1 , provides 
a measure of the field-aligned ion bulk velocity and field-aligned 
ion flux of all ions of a given species which fall within the energy 
passband of the magnetic mass spectrometer. Derivation of the 
field-aligned bulk velocity depends upon the presumption that 
there are three distinct contributors to the species bulk velocity in 
the spacecraft spin plane and frame of reference (sec Figure I). 
These contributors are the spacecraft-induced ram velocity (V r ), 
E x B drift (V KxB ) and field-aligned drift (V||) velocities. The 
spacecraft ram component is obtained from the orbit and attitude 
data base, while the E x B drift component is derived, using meas- 
ured electric and magnetic fields, as 
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where, on the extreme right-hand side. E z is the measured Z 
component of the electric field, expressed in mV rrf and B is 
expressed in gauss. Here, we have used the approximation that the 
magnetic field lies in the spin plane of the DE I spacecraft. As 
was mentioned above, this approximation is always good to 
within 12° for the data used in this study. Then, the total spin 
plane ion species bulk velocity in the spacecraft frame of refer- 
ence is given as 
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where \‘n is the tield-aligned bulk velocity we are trving to deter- 
mine. v rx and i’„. are the X and Y components of the ion ram 
velocity . and 
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Note that positive values ot field-aligned bulk velocity correspond 
to upward and downward drifts in the southern and northern 
hemispheres, respectively. All field-aligned drifts presented in 
this paper arc upward, however, and we have plotted these as 
positive. Here, is parallel to the DE I spin angular 
momentum. The unknown parallel drift velocity is given as 


0| — ( v ry -t- r KxB ) coil a ) - v rx 


(4) 



Spin Phase, <|> 

(degrees) 

Fig. I. The coordinate system and geometry used for computing ion 
species upward bulk velocity from RIMS spin curves and PWI electric field 
data is displayed in the top panel. At the bottom, a typical RIMS count rate 
spin curve is shown, and the meaning of angles identified in the sketch above 
is illustrated. The magnetic pitch angle at which counts arc accumulated is 
displayed as the triangular wave in the bottom panel. 


The angle ot is measured as the difference 

a = <{>,, — <J>b (5) 

in spin phase angles between the location (4>p) of the centroid of 
the distribution of ion counts with spin phase and the direction 
(4 >b) along which the positive magnetic field direction is viewed. 
Spin phase angle (4>) is that measured about +z, from the space- 
craft velocity vector to the RIMS radial head look direction. In 
order to derive the field-aligned ion flux for a given species, the 
product ot the instrument count rate and calibrated instrument 
sensitivity factors are cosine weighted, with respect to the spin 
plane velocity vector (given in equation (2), above), and inte- 
grated over spin phase angles, yielding the flux U T ) along the 
direction of the total spin plane velocity vector. The field-aligned 
flux is then obtained as the product 


h = Jr 


v 'll 
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and the species particle density is estimated as 
n _ J\\ 

l^spin planel ^ 

where |v spm p | alK .| is determined using equation (2). having deter 
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mined Vp vB using measured fields as in equation (I), the ram 
components using orbit and attitude data in conjunction with mau:- 
netic field data, and v j( using equation (4) as described above. Our 
estimate tor the species panicle density is actually a lower bound, 
which disregards the possible presence of ions at energies lower 
than any positive spacecraft potential. This technique has been 
applied to data obtained in and near the 39 upwelling ion event 
intervals included in this study. We have used 1-min data aver- 
ages throughout in analyzing the PWI electric field and the RIMS 
ion data. 

3.3. Error Analysis 

I he above analysis involves two potential sources of error. 
First, tor angles (a) near tt or 0, corresponding to highly field- 
aligned or anti-field-aligned flows, the analysis produces values 
of i'|| which are unrealistically large and positive or negative, 
respectively, for any nonzero difference ( v ^ — |v> xB |). We have 
discarded results which yield a total velocity so large that the 
corresponding ion kinetic energy per charge would be too large to 
be passed by the magnetic mass spectrometer. The accuracy with 
which we can determine a is on the order of ±2°. so that 
the analysis of flows which are in actuality within -2° of being 
tield-aligned would be expected to produce values of v’n whose 
magnitudes are arbitrarily large. We have propagated an uncer- 
tainty of ±2° in ot through the analysis. 

Second, the accuracy with which wc know the values of the Z 
component of the dc electric field is limited by the performance of 
the PWI Z antenna in these low'- to medium-density plasma 
regimes. We have carefully examined the electric field data on a 
case-by-casc basis and selected only those events characterized by 
spin-averaged field profiles which vary smoothly, or with reason- 
ably small amplitude, both immediately prior to and following the 
upwelling event itself, as identified in the ion data. Further, we 
have derived estimates of the magnitude of the uncertainty in the 
values ot E z . in the manner described below, and propagated these 
through the analysis. In processing the electric field data, wc have 
begun by producing 4-s (spacecraft spin period = 6 s> averages 

i U and standard deviations (a 4 J about the means of the 16 
sample s 1 data, producing 15 pairs of values per minute. The 
mean values ot these two quantities represent I -min averages of 
which we have used in the computation of ion bulk 
parameters, and of a 4 s , which represents a spin period time scale 
error estimate. Further, a I -min time scale error estimate (a W) s ) is 
derived as the standard deviation of the 15 values of E z4 s about 
their mean. These two error estimates, in addition to a nominal 
uncertainty ot 4 mV m 1 , have been treated as independent con- 
tributors to the uncertainty in as 

\E Z = Vcri s + crio ,+ (8) 

4. Results 

4.1. Results for a Single Pass 

Figure 2 illustrates the nature of the data used in the analysis 
and the results for a single pass on day 60 of 1982. These results 
are not typical, in that this was a particularly intense event, 
observed at unusually low altitude on a day characterized by a sus- 
tained high level of geomagnetic activity. The ion and electric 
field data shown in Figure 2 arc plotted against UT, with various 
spacecraft orbital parameters also indicated along the horizontal 
axis. The interval selected shows a DE 1 pass from the northern 
polar cap across the day side auroral zone. 


Figure la shows the RIMS radial head spin time spectrogram 
for 0 + and has been described above in the context of data selec- 
tion. The full time resolution Z component electric field data and 
the processed Z component data are shown in Figure lb and 
Figure 2c, respectively. The full time resolution electric field data 
in Figure lb show not only the steady offset induced by spacecraft 
motion across the magnetic field, but also a distinct modulation 
induced by nonstationary spacecraft sheath structure. Recall that 
since the electric Z antenna is aligned with the spacecraft spin 
axis, spin modulation in this signal is not expected in the presence 
of spacecraft sheath structure unless the structure changes with 
spin phase. Note the large-amplitude, unresolved excursions in 
the electric field signature near 2327 UT. The long wire electric 
antenna on DE I shows a similar signature, which is quite 
common near the UWI source region. The possible role of these 
large-amplitude and highly structured electric fields in producing 
the upwelling ion flows will be addressed in a future publication. 
The final electric field data product, which is used in computing 
the ion bulk parameters, is shown in Figure 2c, where values of 
A£, as given in equation (8) are used in determining the size of the 
error bars. In these I -min averages, the contribution induced by 
spacecraft motion has been explicitly removed, and the spin mod- 
ulation removed by the averaging procedure. It is evident from a 
comparison ot Figures lb and 2c that significant structure, present 
in the unprocessed field data, is lost in the I -min data product. 

In Figures Id through If the bulk parameters (V h 7„, and //, 
respectively), derived as described above, are presented for each 
of the ion species H + , He + , CT . and 0 + + . Although the field- 
aligned ion velocities and fluxes are generally upward throughout 
the interval, they are most pronounced during the minute extend- 
ing from 2325:30 UT to 2326:30 UT. It is at this time that the 
event is said to maximize in intensity and that the spacecraft is 
thought to be directly downstream of the upwelling ion flow 
source, in the dayside polar cleft. 

4.2. Typical Flux. Velocity . and Density 
Levels at Event Maxima 

The ramplike upwelling ion signatures observed as the space- 
craft passes from high latitudes through the dayside auroral zone, 
followed by sharp cutoffs at the low-latitude edge, are common 
features ot UWI signatures. These features are interpreted in 
terms of a convective effect, whereby the source of the ion up- 
welling exists in a localized zone near the low-latitude edge of the 
event and ions are convected to higher latitudes, dispersing them 
poleward ot the source. From this point of view, the region where 
the upward flux maximizes, near the low-latitude edge of the 
event, may be considered to be most closely associated with the 
source region. The ions observed in this region are characterized 
by the shortest flight paths from source to observer and, as we 
note, are those with the largest upward field-aligned velocities, 
slower particles having been convected further poleward by the 
time they reach the spacecraft altitude. We have therefore taken 
the ion bulk parameters measured at this location to be indicative 
of the intensity of a given event, with the proviso that there is a 
bias toward larger field-aligned velocities. 

Having selected the analyzed data for 39 events, we have iso- 
lated the ion bulk parameters at the time in each event when the 
maximum flux was observed. Figure 3 shows a scatter plot of the 
measured maximum flux versus the inferred density for each of 
the four species: H ^ . He \ CT . and CT * . Points lying along the 
two lines running through the plot represent upward bulk 
velocities of I km s' 1 and 10 km s 1 . While the measured fluxes 
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2.13 2.06 1.98 1.91 1.84 1.76 1.71 1.65 

9.6 9 9 10.1 10.3 10.5 10.6 10.7 10.9 

19.4 14.8 11.3 8.7 6.9 5.6 4.6 3.9 

69.4 67.0 64.4 61.6 58.6 55.5 52.1 48.7 


Fig. 2. RIMS ion data and PWI electric field data are displayed as functions of time for a dayside auroral zone pass on day 60 of 
1982. The top panel shows a RIMS O’*" spin-time count rate spectrogram in standard RIMS format. The second and third panels from 
the top show PWI dc electric field data in almost raw and processed formats, respectively (see text). The bottom three panels display 
derived ion species bulk parameters of upward velocity, upward flux, and particle density for several ion species. Note that the O ^ ion 
densities have been divided by 10. Various orbital parameters, as well as UT, are listed along the abscissa. 


and densities vary- over several orders of magnitude, the upward 
velocities tend to lie between I and 10 km s _l for the heavier 
species and near or above lOkms 1 for the lighter species. Plots 
of the relative occurrence frequency of maximum density, field- 
aligned velocity, and flux observed in these events for the species 
H , He * , O', and 0 + + are presented in Figures 4a through 
4c, respectively. In the case of the maximum flux (Figure 4c). the 
n ux values plotted along the abscissa have been mapped to a 
common altitude of 1000 km, according to an r } scaling law', 
while the densities and velocities arc plotted as measured, 
irrespective ot spacecraft altitude. This is in contrast to the flux 
\ allies plotted in Figure 3. which are as measured, irrespective of 


altitude. In Figure 4. the species arc differentiated by the plotting 
symbols used. 

Figure 4a shows that these events are typically oxygen-rich 
near the source region with typical fractional and absolute O* 
ion densities of 0.9 and HKK) ions cm respectively. This 
represents somewhat large O* densities for these altitudes 
( — 1.5-2 R F ) as compared with polar cap field lines (see, for 
example, the total electron density results of Persoon ei al. 

| 1 983 1 . or the ion species density results of M. O. Chandler et al. 
(Observations ot polar ion out Hows, submitted to Journal o( 
Geophysical Research. 1989). Such large O ' densities are indica- 
tive ot increased ionospheric scale heights, probably associated 
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Fig. 3. A scatter plot ot measured upward ion flux versus measured ion density at the low-latitude edge of UWls is shown for several 
ion species. The data were gathered over a range of DE I altitudes between 1 .3 and 2. 1 R F . The straight lines running through the plots 
show the loci of points lor ions moving upward at I km s 1 (lower line) and 10 km s 1 (upper line). 



Ion Species Density Species Field Aligned Bulk Velocity 

(km secA-i) 




Fig. 4. Line plots showing the frequency of occurrence of various values 
ot (a) ion density, (b) upward bulk velocity, and (r) upward particle flux, as 
measured at the low-latitude edges of UWIs. The plotted tlux values are 
normalized to a common altitude of 10* km. 
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with bulk ionospheric heating at much lower altitudes. These O ‘ 
density results are consistent with the model polar wind results of 
Si hunk ami Sojka | 1989). Our H * density results are in wide v ari- 
ance with those of Schunk and Sojka. however, with typical 
observed H densities of ^ 100 ions cm ' at geocentric distances 
between 1 .5 and 2 R, as compared with model results in excess of 
several thousand ions per cubic centimeter tut 25(H) km altitude). 

The upwelling 0 + flux (Figure 4< ) is also dominant over that 
associated with the other ions, with maximum scaled flux levels 
which typically approach l() y ions cm : s', far in excess ot 
predicted polar wind 0 + tlux levels. On the other hand, the largest 
field-aligned velocities arc measured in the FI + ions, tor which the 
upward bulk velocity is generally greater than 10 km s _l near the 
latitude of the source region. More modest bulk velocities, ot the 
order of 3 km s" 1 , are typically observed in the ions. 

We have found that the measured upward field-aligned bulk 
velocity varies systematically with species mass per charge, with 
light ions moving upward with higher velocity than the heavier 
ions. This point is illustrated with Figure 5, where characteristic 
upward bulk velocities are plotted as a function of species charge 
per unit mass. The distribution of velocities shown in Figure 5 is 
dominated by the FI ~ ions (at qim = I) which arc clearly moving 
upward much faster than the other species, although the He" ions 
are also significantly faster than the two species (0 + and O" ") 
plotted. The trend is not as strongly supported when comparing 
the O ’ and O" + ions either in this figure, which is derived from 
the mean values of the bulk velocities observed at maximum tlux, 
or on a casc-by-casc basis, where the 0 + is often seen to be 
moving upward with higher velocity than the + . We have 
shown a power law r fit to the data in Figure 5, which yields 

~ n \ Z (°‘ 6 km 

1 11 | (mlm p ) | s 

where Z is the ion charge state and (m/m p ) is the ion mass, norma- 
lized to the proton mass. Processes yielding energies proportional 
to ion charge hut independent of ion mass would give rise to a 
square root dependence of v’n upon Uf/m). An example of such a 
process is acceleration in a parallel (to B) electric field. However, 
the warm and anisotropic nature of these upwelling plasmas, 
along with their high density, suggests that the outflow arises due 
to transverse heating in the topside of an already warm cleft iono- 


Fig. 5. Mean values of ion species upward bulk velocities are plotted 
versus the species charge per unit mass. The velocities plotted are those 
derived from measurements taken at the lower-latiiude edges of UWIs, 
where the upward flux is seen to maximize. The line indicates a power law t it 
to the daia. with the fit parameters indicated. 


sphere, followed by upward magnetic folding of the ion dis- 
tributions. The analysis of the thermal properties of these events is 
beyond the scope of this report and will be reported in a future 
paper. However, it is to be noted that upwelling ion flows are 
typically warm and anisotropic, with temperatures (transverse to 
B) of several electron volts, as reported by Moore et ai. [ I986/>|. 

Defcourt ct a/. | 1989| have modeled the evolution of these ion 
flows using particle ray tracing techniques in realistic gravita- 
tional. electric, and magnetic fields. They have found that the up- 
welling ions generally escape the Earth's immediate vicinity in 
large numbers and travel in an extended plume over the polar cap 
into the geomagnetic tail and plasma sheet, where the heavier 
species are rather dramatically accelerated as they cross the 
neutral sheet. 

4.3. Bulk Parameter Correlations and Dependence 
on Location and Geomagnetic Activity 

The relationship between spatial location and the occurrence of 
ion upwelling has been investigated by Lockwood et al. 1 19851. 
Their study showed that for DE I altitudes below 3 /?*.> geocen- 
tric, the probability of observing the asymmetric RIMS spin-time 
signature of an upwelling ion event peaks at near 0.7 in the 
prenoon magnetic local time (MLT) sector and at invariant lati- 
tudes < 1L) between 75° and 80° (see their Plate 4c). The question 
arises as to whether the intensity of the upward flow shows 
systematic variation with the location of the source region. In 
Figure 6, we present scatter plots of the upwelling O" tlux, meas- 
ured during the minute of maximum intensity within the event and 
normalized to a common altitude of 1000 km, versus MET 
(Figure 6a) and IL (Figure 6/;). There is substantial scatter in the 
observations, with several particularly large flux values, and no 
apparent systematic trend. The intensity of the outflow' does not 
seem to be correlated with the latitude or MLT at which the obser- 
vations are made within the source region. 

There is. however, some indication that the ion composition 
within these events varies systematically with magnetic local 
time, with the fractional minor species density increasing toward 
local noon. Figure 7 provides a convenient summary of the ion 
composition for a number of cases in which all four ion species 
were measured simultaneously. As noted earlier, the bulk ot the 
ion plasma consists of O". comprising about 9(Y7r of the ion 
density, while H" contributes — 1 0*55- . and He + and ^ each 
contribute ~l c /c. The upwelling flux typically represents a 
smaller percentage < ■ — 70%) of the total upward tlux, due primari- 
ly to the larger H'*’ bulk velocities. In Figure 7. we show the frac- 
tional ion composition plotted versus the magnitude iMLT— 12| to 
illustrate the existence of a possible systematic relationship 
between proximity to local noon and ion composition within these 
events. Note that values are plotted along the abscissa without 
regard to whether data were collected in the prenoon or postnoon 
sector. The fractional composition of the minor species appears to 
increase with proximity to local magnetic noon. The case for the 
compositional dependence upon MLT is stronger for O " * (linear 
correlation coefficient = -0.7) than for the other species. The 
reason for this dependence is unclear. The O" density (and 
therefore total ion density) appears to decrease with proximity to 
noon, although the data arc scattered (linear correlation coeffi- 
cient |tf|— 0.3), while the absolute minor species densities show 
no significant trend. These apparent trends in the ion composition 
could be anomalous, resulting from correlations between MLT 
and altitude imposed by the DE 1 orbit. Further study of this 
subject will be required. 
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As shown by Lockwood el at. |I985), the invariant latitude c 
(he upwelling ton source region varies systematically wit 
geomagnetic activity, moving to lower latitudes with inercasin 
values ot the 3-hour Kp index. The association of this upwellin 
source region with the geomagnetic polar cleft and its systemati 
latitudinal variation is further illustrated in Figure 8. where th 
invariant latitude at which the maximum upwellins O + flux occur 
is plotted versus a measure of the interplanetary maenctic ficl. 
(IMF) B z component. The values of B. plotted were obtained fron 
hourly averages published by the National Space Science Dan 
Center (NSSDC) World Data Center A for Rockets and Satellites 
l hese data are a compilation based on IMP 8 and ISEE 3 observa 
turns. The plotted values represent weighted (and corrected foi 
solar wind transit time) hourly averages of B. during the houi 


immediately preceding the minute in UT during which the flux 
observation was obtained. For example, in the case of day 71 of 
1982 the maximum O* flux was observed at 1934 hours UT. The 
value ot B. associated with that event is taken as 

rxim-isKMol 

' " 60 < 9 » 

where was obtained from the World Data Center data 

base. We have indicated a least squares linear fit to the data, 
which yields a correlation coefficient of 0.72. We have also 
shown (dashed lines) the linear least squares results of Carbary 
and Meng I '986|. giving the fitted locations of the poleward and 
equatorward cusp boundaries, based on DMSP F-2 and F-4 low- 
energy electron signatures, as functions of B : (measured by the 
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Fig. 8. A scalier plot of the latitude at which the upward CU flux maxi- 
mizes within UWIs. versus the preceding hourly average of the IMF B 
component (corrected lor solar wind transit time to the magnetopause). The 
solid line represents a linear least squares fit to the data. The two dashed I ines 
show the locations of the upper- and lower-latitude edges of the cleft! ike (see 
text) electron precipitation signatures as functions of B.. These edges are as 
determined in the study of Carbary atui Meng [1986], 


IMP 8 /GSFC (Goddard Space Flight Center) magnetometer). We 
note here that although Carbary and Meng ( 1986] use the word 
“cusp." the study they performed applies to a broader region (in 
MLT) which is more properly referred to as the "magnetospheric 
cleft" f Carbary and Meng , 1986. Figure 2 ; Newell and Meng , 
1 988 1 . The distinction is significant because the cusp, found 
within a region highly localized near noon MLT, as thought to 
connect to the dayside magnetopause, while the cleft, found 
within a broad region which stretches across the dayside, is 
associated with the low-latitude boundary layer. Although, as in 
the Carbary and Meng study, there is significant scatter in the 
data. Figure 8 unambiguously locates the upwelling ion source 
region within the magnetospheric polar cleft, as identified by the 
low-energy electron signature. 

There does not, however, appear to exist a strong relationship 
be tween cither the occurrence probability or intensity of up- 
welling ions and the averaged value of the IMF B z during the hour 
immediately preceding the observation of ion upwelling in the 
RIMS data. In Figure 9 . we study the occurrence of values of B , as 
observed, both in conjunction with an UWI observation and 
regardless of UWI observation. All three panels in Figure 9 
feature an abscissa which is composed of B z values binned in I y 
bins, trom - 6 y to +6 y, with the leftmost bin including all 
observations of B z less than -6 y and the rightmost bin those with 
B z greater than or equal to + 67 . There were a total of 70 
observed UWI events for which the IMF data were available dur- 
ing both the present and previous hours, allowing the weighted 
average described in the previous paragraph to be performed. The 
value plotted along the ordinate ot Figure 9 a represents the 
percentage ot those 70 events which were characterized by a 
weighted hourly average of B z which fell within the given bin. 
The error bar is derived from the square root of the number of 
samples in the bin. For example, 9±3 of the events were asso- 
ciated w ith a value of B z of less than —6 7 , yielding 12.8 ± 4 . 3 % 
ot the 70 samples. In Figure 9b the occurrence of frequency of B, 
values, tor all available hourly IMF data during the same time 
Period as in Figure 9 a (days 36^268 of 1982), are plotted for 
comparison. The distribution in Figure 9b is peaked near B . = 0 


flf 
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Fig. 9. The occurrence, during the period between days 36 and 268 of 
1982, of various values of the IMF B z component is studied, both in relation- 
ship with the occurrence of UWIs and independent of such occurrence. <«) 
The distribution of B. values measured during the hour immediately preced- 
ing a RIMS UWI observation, {b) The distribution of A. for all the available 
data, without regard to RIMS observations, (e) The difference between the 
data in Figures 9 a and 9b. These differences should reveal any significant 
trend in the dependence of the occurrence of dayside ion upwelling upon B . 
All IMF data are corrected for solar wind transit time to the magnetopause 


and falls off quite symmetrically to either side. The large 
occurrences in the extreme bins are a reflection of the larger bin 
size. Figure 9c displays the differences between the data in the 
above two panels, such that the plotted values are equal to those 
plotted in Figure 9a, less those plotted in Figure 9b. If the occur- 
rence of UWI events were dependent upon B z . we would expect to 
see that reflected most clearly in Figure 9c. If, for example. UWIs 
were more likely to occur for negative values of B z . we would 
expect to see positive differences in Figure 9c for B z < 0 and neg- 
ative differences for B z > 0 . No such trend is evident, and we 
conclude from these data that UWI occurrence is not dependent 
upon the average value ot B z during the hour immediately preced- 
ing event observation. 

A similar statement can be made regarding the dependence of 
upwelling 0 + intensity on the averaged B. values, as illustrated 
using the scatter plot of Figure 10 . Here we have plotted the 
maximum upwelling CT flux, normalized to 1000 km altitude, 
measured during a given event versus the hourly averaged value 
of B z ^r the hour immediately preceding event observation. The 
linear correlation coefficient of -0.14 indicates no significant 
dependence among the parameters plotted. These results are 
somewhat surprising, in view of recent observations of Lockwood 
et al. (1988] which suggest that dayside ion upwelling may be a 
direct consequence of magnetic coupling at the subsolar magneto- 
pause. Further study of the relationship of UWI events to 
interplanetary and magnetospheric conditions will be the subject 
of a future publication. 

4.4. Cleft Region Source Strength 

The upward ion flux levels discussed in the previous sections 
arc not typical of the entire dayside auroral zone, but are 
representative of the most intense portions of the upwelling flows. 
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Fig. 10. Upwelling O" flux is displayed as a scatter plot versus the 
immediately preceding hourly averaged value of the IMF B. component, 
corrected for solar wind transit time to the magnetopause. The flux values 
plotted are the maximum values measured within the UWIs and are normal- 
ized to a common altitude of I0 3 km. A linear fit to the data, as well as the 
resulting linear correlation coefficient, is indicated. 


in the immediate vicinity of the source region. We have computed 
bulk parameters from regions downstream (generally poleward) of 
the source region as well. All the measured fluxes, as well as a 
measure of the occurrence frequency of davsidc ion upwelling. 
need to be incorporated into any estimate of the total ion outflow 
from this ionospheric source. The expectation value for the up- 


ward particle flux at any location may be expressed as the product 
of the upward flux measured at that location at the time when an 
upwelling ion event is observed and the probability that an up- 
welling ion event will be observed at that location. We have 
binned our ion flux measurements in IL and MLT and combined 
them as a product with similarly binned values of the occurrence 
probabilities reported by Lockw ood cr a ! . 1 1985|. We have used 5° 
IL bins between 5() c and 90' and 1-hour MLT bins from 0 to 24 
hours. The mean flux in a given bin is taken as 

N 


where J ij n represents the mh of a total of N flux observations 
(normalized to 1000 km altitude) in the angular bin associated 
with the indices (/./). The expectation value for the flux in the bin 
(i.j) is then given as 

~ (A ij) x if ij) (11) 

where / is the occurrence probability reported by Lockwood ei 
cd. 1 1985], Expectation values computed in this manner are 
presented in Figure I I for the species FP , Fie ~ . O" . and " . 
where the value of logj,,t<y>,-j) is plotted in polar grey scale 
versus IL and MLT. The expectation values plotted in Figure 1 1 
represent the closest estimate we can make at this time of the 
mean upwelling source strength for the various ion species 
shown. The outflowing source region may extend poleward of the 
region shown in Figure 11: however, such an extension is not 
identifiable within the context of the UW1 signature identification 
on which this work is based. As in the case of the maximum event 
fluxes reported in previous sections. O* and FI* show the most 
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Fig. 11. The expectation value <7 !( > for the upward parallel flux is plotted for the ion species H ‘ . O * , O + ^ . and He * . in polar 
grey scale format versus IL and MLT. The flux in a given 1L-MLT bin is derived as the product of the average flux measured in that bin 
within UWIs and the probability that an UWT will be observed within the bin. 
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intense outflows, followed by He* and CT T In ail eases, the 
source region is rather spatially structureless. 

The fluxes shown in Figure 1 I represent the ion flux expected 
at an altitude of 1000 km. having been sealed from the observa- 
tion altitudes according to an (r/r n ) 3 dependence, where r and 
arc the geocentric radii of the observation point and the 1000 km 
altitude surface, respectively. The total hemispheric particle out- 
flow rate may be obtained by integrating over the distributions 
shown in Figure 11. To facilitate the integration, the fluxes are 
further scaled to the Earth’s surface, where IL and the spherical 
magnetic latitude are equivalent. The total terrestrial outflow (/) 
from this dayside source region, for a given ion species, is then 
given as 

1 = 2 f <J(/\\,<J>)> da 
n 

y r2v nix 

jf) <*<(> dk <Ji{r 0 ,k.4>)> ( 1 2) 

where we have approximated the geomagnetic field lines as radial 
and the factor of 2 accounts for source regions in each 
hemisphere. This integration has been numerically carried out 
over the distributions of fluxes shown in Figure 1 1 . The resulting 
outflows are listed in Table 1. These numbers may underestimate 
the total outflow, since significant UW1 outflows may occur pole- 
ward of regions identified as being within UWIs, 


TABLE I. The Total Rate <J\\> at Which Ions Flow 
out of the Dayside Cleft Region 


Ion Species 

<7(i>. ions s 1 


5 x |(F J 

He^ 

9 x |(F 

cr 

2 x 10- 

cr 

2 x I0 2S 


Flows from both hemispheres are included. 


The ion outflow is dominated by CT , with 2 x I0~ ions s 1 
emanating from the dayside cleft region. The H* source strength 
represents —20 c /c of the total outflow, while the minor species 
{O* * and HeM arc at the several percent level or less. These 
results are to be compared with those of You et al. 1 19851 who 
studied major ion outflows at energies such that 10 V ^ E/q ^ 17 
kV. They reported extensively on outflows at latitudes including 
the auroral zone and polar cap. distinguishing among flow levels 
on the basis of geomagnetic activity, MLT/IL regime, and solar 
cycle. Many of the results of You et id. 1 1985] are summarized in 
their Table 1 . They found the total net auroral zone ion outflow to 
lie between 2 x I0 -5 and 4 x I0 25 ionss" 1 for and between I 
x 1 O' 7 ' and 8 x 10 2s for CT , depending on geomagnetic activity 
and solar cycle. In our Table I , we show a total UW1 O* outflow 
of 2 x 10- ions s' 1 . This is of the same order as that reported by 
Yau et al. during geomagnetically quiet times. On the other hand, 
our results show a total outflow of 5 x 10 24 ions s _l , a factor 
of 5 smaller than that reported by Yau et al. during gcomag- 
netically quiet times. We believe these figures apply to ions with 
energies less than or of the order of 10 c\lq. This is based on both 
the published work of Moore et id. [1986/?| and on other, 
unpublished work, which show typical ion thermal energies 
within these events of 5-10 cV/q. The flux levels in the current 
study arc representative of a wide range of geomagnetic activity. 
The available data do not permit breaking down the global UWI 



outflow according to geomagnetic activity. It appears, especially 
in the case of the 0* . that the UWI source represents a substan- 
tial fraction of the total global ion outflow. 

5. Conclusions 

With the data presented in the sections above we have reported 
the bulk properties of the upwclling ion plasma in the dayside 
cleft region. We have quantified the distribution of ion species 
density, upward bulk velocity, and flux within the cleft source 
region. Additionally, we have combined the bulk flow' results 
from a number of these events with the occurrence probabilities 
obtained by Lockwood et al. 1 1985] to give reliable estimates of 
the source strength distribution as a function of MLT and 1L. 
Finally, we have presented a measure of the total cleft region 
source strength in terms of the number of ions per second to flow 
away from the ionosphere. The following specific conclusions 
may be drawn from the work presented above. 

1. Ail observed ion species are seen to participate in the ion 
upwclling. Within a given event, the upward flux typically max- 
imizes near the low-latitude edge of the outflow. At that location, 
typical ion species densities, upward field-aligned bulk velocities, 
and upward field-aligned fluxes (mapped to a common altitude of 
10(X) km) are summarized in Table 2. The distributions about 
these typical values have been presented in Figure 4. 


TABLE 2. Typical Ion Bulk Parameters as Measured Wiihin. but Near the 
Lower-Latitude Edge of Upwclling Ion Events 


Species Bulk Parameter 

O' 

H * 

He* 

cr - 

/?. cm ■' 

600 

50 

15 

7 

ru, km s 1 

3 

16 

7 

4 

7 f |,* particles cm 7 s 1 

6 x |0 K 

2 x 10* 

l() 7 

6 x I0 (1 


^Scaled to a common altitude of I (XX) km. 


2. There is a systematic dependence of the upward field- 
aligned ion bulk velocity upon the species charge-to-mass ratio, 
with the velocity increasing approximately with the square root of 
q m. This dependence is illustrated in Figure 5. 

3. The intensity of UWIs, as measured by the maximum up- 
ward CT flux observed within an event, does not show significant 
systematic variation as a function of MLT or IL. This is indicative 
of a rather uniform outflow from the dayside cleft region. This 
point is illustrated in Figure 6, where scatter plots of the max- 
imum normalized CT flux versus MLT and IL arc presented. The 
point is further illustrated, with respect to all measured flux values 
for the species H * , Hc^, (T, and CT in Figure 11, where 
the expectation value of the normalized upward flux is shown 
plotted in polar form. 

4. The ion density composition of the outflow shows a weak 
dependence upon MLT. however, with minor species, particular- 
ly + , representing a larger fractional composition near magnetic 
local noon. This point is illustrated in Figure 7. 

5. The statistical association of UWIs with the dayside mag- 
netospheric cleft, as identified by low-energy electron signatures 
[Carbary and Men% % 1986], is illustrated in-Figure 8. This associ- 
ation points up the likelihood of a relationship between the 
phenomenon of dayside ion upwelling and the solar wind- 
magnetosphere interaction. 

6. There is, however, no clear-cut relationship between either 
the occurrence or intensity of UWIs and the immediately preced- 
ing hourly averaged value of the IMF B z component. These facts 
are illustrated in Figures 9 and 10, respectively. As the simplest 
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mea s ur c ot the degree of electromagnetic coupling between the 
magnetosphere and solar wind, one might expect to find increased 
occurrence or intensity ot UWIs with increasingly negative B. if, 
tor example. UWIs were directly associated with variable aspects 
ot the large-scale auroral convection morphology or with the 
occurrence ot magnetopause flux transfer events, both of which 
depend upon B z . The results presented here indicate that neither is 
the case. 

7. Our estimate of the ion outflow rates (including both 
hemispheres) trom this dayside cleft source have been presented 
in Table I. These estimates arc given based on data collected 
between days 36 and 137 of 1982, and almost all were from 
northern hemisphere passes. Thus, they are generally representa- 
tive of spring conditions during an active portion of the solar 
cycle. Virtually all levels of geomagnetic activity are represented. 
In the case ot 0 + and H^, these numbers validate the cleft ion 
tountain contributions used by Chappell et al. [1987] and 
Delcourt et al. [ 1 989] in demonstrating the significance of the 
ionosphere as a magnetospheric plasma source. 
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Statistical survey of pitch angle distributions in core (0- 
50 eV) ions from Dynamics Explorer 1: Outflow in the 
auroral zone, polar cap, and cusp 

B. L. Giles, 1 C. R. Chappell, 2 T. E. Moore, 1 R. H. Comfort, 3 and J. H. Waite Jr. 4 

Abstract. Core (0-50 eV) ion pitch angle measurements from the retarding ion mass spec- 
trometer on Dynamics Explorer 1 are examined with respect to magnetic disturbance, invariant 
latitude, magnetic local time, and altitude for ions H + , He*, 0\ M/Z=2 (D* or He**), and O** 
Included are outflow events in the auroral zone, polar cap, and cusp, separated into altitude re- 
gions below and above 3 R E . In addition to the customary division into beam, conic, and up- 
welling distributions, the high-latitude observations fall into three categories corresponding to 
ion bulk speeds that are (I) less than, (2) comparable to, or (3) faster than that of the spacecraft. 
This separation, along with the altitude partition, serves to identify conditions under which iono- 
spheric source ions are gravitationally bound and when they are more energetic and able to escape 
to the outer magnetosphere. Features of the cleft ion fountain inferred from single event studies 
are deariy identifiable in the statistical results. In addition, it is found that the dayside pre-noon 
cleft is a consistent source of escape velocity low-energy ions regardless of species or activity 
level and the dayside afternoon cleft, or auroral zone, becomes an additional source for increased 
activity. The auroral oval as a whole appears to be a steady source of escape velocity H*, a 
steady source of escape velocity He* ions for the dusk sector, and a source of escape velocity 
heavy ions for dusk local times primarily during increased activity. The polar cap above the au- 
roral zone is a consistent source of low-energy ions, although only the lighter mass particles ap- 
pear to have sufficient velocity, on average, to escape to higher altitudes. The observations sup- 
port two concepts for outflow: (1 ) The cleft ion fountain consists of ionospheric plasma of 1-20 
eV energy streaming upward into the magnetosphere where high-latitude convection electric 
fields cause poleward dispersion. (2) The auroral ion fountain involves field-aligned beams 
which flow out along auroral latitude field lines; and, in addition, for late afternoon local times, 
they experience additional acceleration such that the ion energy distribution tends to exceed the 
detection range of the instrument (>50-60 eV). 


1. Introduction 

At polar latitudes, because the Earth’s dipole field is nearly 
vertical, charged particles in the ionosphere created by pho- 
toionization or particle impact ionization are relatively free to 
move upward along the magnetic field lines and populate 
higher-altitude regions of the magnetosphere. Since the dis- 
covery of energetic CP ions by Shelley et al. [1972], ion com- 
position measurements have revealed that a significant 
amount of plasma in the Earth’s magnetosphere is of iono- 
spheric origin [e.g., Chappell et al ., 1987; Moore et al, 
1989; Moore , 1991]. The importance of this terrestrial source 
is confirmed by the presence of energized ionospheric ions 
throughout the magnetosphere, such as the ring current and 
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plasma sheet [Johnson et al 1977; Geiss et al ., 1978; 
Lennartsson et al, 1979; Balsiger et al ., 1980; Lundin et ail 
1980, Peterson et al., 1981; Sharp et al. , 1982; Eastman et 
al., 1984, Lennartsson and Shelley, 1986], the magnetopause 
and magnetosheath region [Peterson et al., 1982; Hultqvist , 
1983], the magnetotail lobes [Sharp et al . , 1981], and the 
plasma mantle and magnetotail boundary layer [Frank et al., 
1977; Hardy et al . , 1977; Lundin et al., 1982a, b; Candidi et 
al., 1982]. In addition, it has been shown that a significant 
portion of these ionospheric-source ions have energies below 
1 keV [Yau et al., 1934; Giles et al, 1988; Delcourt et al, 
1988]. These lower-energy particles are very sensitive to 
small fluctuations in local magnetic and electric fields. They 
interact with charge carriers for magnetospheric and iono- 
spheric currents and operate as a permeating medium in which 
plasma waves are generated and propagated to other parts of 
the magnetosphere thereby influencing larger-scale interac- 
tions. Thus information about the dynamics of this plasma 
component is an effective tracer for broader scale magneto- 
spheric processes. 

In understanding the transport and dynamics of charged par- 
ticles, information about the particle pitch angle distribution 
(PAD), the arrangement of directional flux with respect the lo- 
cal magnetic field direction, is an invaluable aid. For exam- 
ple, case studies by Sharp et al. [1977, 1979] exploited field- 
aligned pitch angle distributions to infer the presence, and 
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vertical extent of, parallel potential drops within an auroral 
acceleration region. Also, occurrence statistics of distinct 
PAD classifications have been used to examine possible 
plasma source locations, acceleration regions, and the subse- 
quent evolution of plasma distributions during transport [e.g.. 
Nagai et al, 1983; You el al.. 1984], The lowest-energy ions,’ 
in particular, have been shown to exhibit complex pitch angle 
characteristics such as field-aligned beams, conics, and 
trapped (pancake) distributions [Chappell, 1982; Chappell et 
al., 1982a], 

The retarding ion mass spectrometer (RIMS) experiment on 
Dynamics Explorer I (DE 1) represents a major effort to thor- 
oughly characterize the lowest-energy plasma component by 
measuring ions ranging from 0 eV (or spacecraft potential) up 
to 50 eV throughout almost the entire volume of the inner 
magnetosphere, from the topside ionosphere out to 4.67 R E . 
The RIMS effort has, in its lifetime, revealed previously unob- 
served ionospheric plasma sources to the magnetosphere and 
called into question accepted theories concerning the relative 
contributions of the internal and external sources to magneto- 
spheric plasma content. Specific to the RIMS instrument are 
observations of the supersonic polar wind [Nagai el at., 
1984], a source of low-energy heavy ions from the polar cap 
and the cleft region [Waite el aC, 1985; Lockwood et al., 
1985a, b], upward fluxes of N*. N** [Chappell et al.. 1982b] 
and of molecular ions N 2 *, NO*, and 0 2 * [Craven et al., 1985], 

It is this data set that will be used to obtain pitch angie distri- 
butions ot the individual charged particle species. 

One of the most intriguing topics in magnetospheric 
physics research that will be addressed is the question 
regarding the source of particles populating the inner magne- 
tosphere. Concentrations and relative dynamic behaviors of 
H , He**, and O* are commonly used as indicators of a solar 
wind or ionospheric source. The RIMS mass spectrometer is 
capable of resolving several mass peaks, including H*, He*, 
O*. O**, and one at a charge-to-mass ratio of 2. Because RIMS 
provides no direct information on charge state, these latter 
measurements could conceivably be H 2 *, He**, or D*. Breig et 
al. [1987] and Breig and Hanson [1991] have reported on deu- 
terium ion flows in the thermosphere. In this case the instru- 
ment is a predecessor of RIMS which again, yields the M/Z=2 
mass peak with no charge information, but additional correla- 
tive evidence with H* determines these samples to be D*. 
Young et at. [1977, 1982] and Comfort et al. [1988] report 
observations in the relatively cold outer plasmasphere that are 
most likely terrestrial source He**. In the present work, be- 
cause the DE 1 orbit samples low altitudes where D* measure- 
ments are likely and samples higher altitudes where low- 
energy He** measurements are possible, the RIMS 
measurements will be referred to as being of the M/ Z=2 ion. It 
is thought, because of the RIMS energy range, that any He** 
measurements would represent doubly-ionized helium produced 
in the ionosphere or plasmasphere, not cooled solar wind 
He**. Regardless though, of whether the Af/Z=2 measurements 
are representative of He** or D*, the key issue is whether 
evidence exists for these ions, and of H*, He*, O*, and O**, to 
have an ionospheric source and then flow outward, such that 
they populate both inner and outer portions of the 
magnetosphere. 

In section 2 we describe the RIMS instrument and the DE 
orbit, and outline the pattern recognition technique used. 
Section 3 describes the classification of the data in terms of 
pitch angle distributions and details the consequences of ob- 
serving low-energy particles with a moving spacecraft. 


Section 4 contains occurrence probability data with respect to 
magnetic disturbance, invariant latitude, magnetic local time 
(MLT), and geocentric altitude for invariant latitudes above 
60°. Section 5 discusses these data in relation to the transport 
and dynamics of plasma in the Earth's inner magnetosphere. 
Section 6 is a summary. A paper (Giles et al., Statistical sur- 
vey of pitch angle distributions in core (0-50 eV) ions from 
Dynamics Explorer 1: Map of the magnetosphere’s low- 
energy ion component, unpublished manuscript, 1994) 
presents statistics for invariant latitudes below 60°, 
examining trends with disturbance, magnetic latitude, MLT 
and L shell. 

2. Instrumentation and Data Analysis 
Considerations 

The survey covers the ions routinely measured by the DE 
1/RIMS instrument (H + , He + , 0 + , Af/Z=2, and 0 + +) in the energy 
range 0-50 eV. The observed count rates for each ion were 
classified according to their pitch angie distribution as deter- 
mined Irom the radial sensor head of the instrument as it ro- 
tated through a full range ot pitch angles in each 6-s satellite 
spin period. As the spacecraft spins, the radial head (mounted 
perpendicular to the spin axis) responds to ion fluxes for 
nearly the full range of angles (-180° to +180°) with respect to 
the spacecraft ram direction. With knowledge of the Earth's 
magnetic field in relation to spacecraft ram, the RIMS spin 
angle response corresponds closely to a pitch angle distribu- 
tion (detailed in section 3). The instrument and its operational 
characteristics are described by Chappell et al [1981], Fields 
et al. [1982], and Olsen et al [1985]. 

The elliptical polar orbit has perigee of 675-km altitude and 
apogee of 24,875-km altitude. The motion of the Earth around 
the Sun causes the orbit plane to drift westward at a rate of 1 
hour MLT every 15.4 days (covering 24 hours in 12 months); 
the orbit line of apsides drifts 0.328° geographic latitude each 
day (moving from one pole to the other in 18 months). AH 
available data, from October 19, 1981, to December 31, 1984 
(a full precession/rotation cycle), are used to provide nearly 
complete local time and latitudinal coverage out to 4.67 R E . 
However, the 3 to 1 ratio between the drift of the line of ap- 
sides and the westward drift of the orbit plane does result in 
specific areas of poor sampling in local time for a given alti- 
tude. For example, in the 0 to 6 MLT sector, midaititudes in 
the auroral latitude regions are not well sampled. In addition 
to regions not reached by the orbit, the instrument is not al- 
ways operated down to the DE 1 perigee of 675 km. At low al- 
titudes the plasma density can be quite high and the instrument 
often enters a shutoff mode to protect the particle-detecting 
channeltrons from excessive count rates. 

The instrument operational modes allowed programmed 
sampling of the mass range such that the major masses were 
sampled more frequently and constitute a larger portion of the 
database. Table 1 lists the number of 1-min spin angle distri- 
bution samples for each species above 60° invariant latitude 
(ILAT) in four MLT quadrants and two geocentric distance lev- 
els. Over 500,000 samples were obtained for the H\ He + , and 
° + ions and over 400,000 samples for the M/Z=2 and 0 ++ ions. 
Because of the varying sample rate, the statistical calculations 
are carried out for each ion species separately; the probabili- 
ties presented are characteristic of the population of individual 
ion species, not of the ion population as a whole. Fifty-eight 
percent ot the data were obtained during periods of relatively 
low activity ( Kp < 3) and 42% during high activity ( Kp > 3). 
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It should be noted that, for the 0 to 50 eV energy range 
floating spacecraft potentials can prevent the detection of the 
lowest-energy ions. The RIMS team dealt with this spacecraft 
charging by the addition of a biasable voltage plate at the en- 
trance aperture programmable to 0, -2 -4, or -8 V relative to 
the spacecraft potential. The present data set includes all 
available data, regardless of aperture bias setting. It is possi- 
ble that the data set does not include particles with ram energy- 
per-charge less than the potential of the spacecraft 

For each ion the raw RIMS data are read and the counts-per- 
accumulation period (12 ms) are averaged over the time resolu- 
tion period for each spin angle bin. For this study, I -min time 
averages with 5° spin angle bins were chosen. This 5° spin 
angle bin refers to the center of the view cone and does not 
imply 5 resolution in the measurements. A subroutine scans 
the counts-versus-spin-phase-angle array to identify and char- 
acterize, in terms of centroids, widths, and skews, local max- 
ima m the distribution with respect to the magnetic field direc- 
tion. Next, a filtering pattern recognition routine produces 
angle distribution description for each minute of 
RIMS data. When combined with orbit parameters and mag- 
netic activity indices, this is the source for statistical analy- 
sis. Further description of the data analysis may be found in 
the work by Giles [1993], 

The relations used to calculate the centroid and coefficient 
of skewness, taken from Beyer [1984], apply to discrete data 
grouped into a distribution having spin angle bin classifica- 
tions 0,, i _ /, 2, .... k, and corresponding measurement val- 
C3Se ’ instrum ent count rate per spin angle bin), 
t- 1.2 k. The centroid is an arithmetic mean given 


where A/ 2 and are the second and third moments about the 

mean of the sample. The rth moment about the mean, M is 
given by 

k 

sf n* \f n* ^ \ r 
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Da,a in Terms of 

3.1. Relationship to Pitch Angle Distributions 

Strictly speaking, spacecraft spin angle distributions ob- 
served m situ are no. true pitch angle distributions as cus.om- 
nly defined. The measured spin angle is defined in the refer- 
ence frame of the spacecraft, which is different from the 
plasma rest frame due to spacecraft motion (2 to 8 km s 1 de 
pending on spacecraft altitude). This difference in reference 
l * tU * conse£ i uence Provided the thermal speed. 
( Jit) , of the plasma is large compared to the spacecraft 
ram speed. For an average spacecraft speed of 5 km s ' a fac- 
tor-of-5 ratio of thermal speed to spacecraft speed requires 3 5 
eV for FT, 1 3 eV for He + , and 52 eV for 0 + . 

Because RIMS responds to ions with energies considerably 
below this threshold, the RIMS spin angle data contain infor- 
mation about both the angular distribution of the plasma in its 
own frame of reference (the conventional pitch angle distribu- 
tion) and the drift motion of the plasma relative to the space- 
craft. These two aspects of instrument response tend to domi- 
nate independently - one at times obscuring the other - de- 
pending on the plasma's thermal speed and its relation to the 
spacecraft and plasma drift speed. Figure 1 illustrates the ex- 
tremes of this motional influence. 

Figure I shows simulated RIMS spin angle responses to a 
bi-Maxwellian plasma distribution (T ± /T H = 5) for four combi- 
nations of plasma thermal and drift speeds. The simulation 
adopts a numerical quadrature approach to the solution of the 
response integral, using a sum over aperture area, solid angle 
and energy and weighted with the appropriate instrument re- 
sponse factor to obtain the correct energy-angle dependent re- 
lationship between detector count rate and directional flux 
The parameters illustrating the limits of response were chosen 
or ease of understanding but without strict adherence to phys- 
ical likehood in any given region of the polar ionosphere In 
each case the spacecraft speed is set at 5 km s ' perpendicular 
to B. Shown beside each spin angle plot are two phase-space 
schematic diagrams to clarify the orientation of the various 
spacecraft and plasma components. Both display velocity 
space, one in the spacecraft reference frame, the other in the 
frame of reference of the plasma. 

In Figure la the H* ions have the perpendicular thermal 
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Figure 1. Simulated retarding ion mass spectrometer (RIMS) spin angle responses to bi-Maxwellian 
plasma distributions (T’j^/7'jj = 5). Also shown are phase space diagrams, drawn schematically, locating the 
various spacecraft and plasma components. Note that instrument response is to ions with velocities anti- 
parallel to the detector look direction, (a) H + ions with perpendicular thermal speed component greater than 
spacecraft speed by a factor of 5. The top panel assumes a stationary plasma and the bottom panel assumes 
Vo = 25 km s'\ parallel to B. (b) 0 + ions with perpendicular thermal speed a factor of 5 below the spacecraft 
speed and with same drift speed assumptions as Figure la. For plasmas with thermal speeds smaller than the 
spacecraft and drift speeds, relative drift motion dominates the response. For plasmas with thermal speeds 
larger than the spacecraft and drift speeds, the true plasma frame distribution function emerges. 


speed component greater than the spacecraft speed by a factor 
of 5 (V th - 25 km s* 1 and V sc = 5 km s' 1 ). In the top panel the 
RIMS response is to a stationary plasma (drift speed V D = 0); 
the angular characteristics of the distribution are clearly visi- 
ble as peaks in the spin angle directions perpendicular to the 
magnetic field (0° and 180°). The relative motion of the 
spacecraft (5 km s' 1 with respect to the stationary plasma) ap- 
pears as an elevation in count rate in the direction of space- 
craft motion (0°). The bottom panel assumes V D = 25 km s' 1 , 
parallel to B. In this case, bulk drift speed comparable to 
thermal speed, the angular characteristics of the distribution 
have become nearly unresolvable although the additional drift 
speed of the distribution may be identified by the shift from 


ram in the direction corresponding to an addition of the space- 
craft and drift velocity vectors. (Without regard to the rela- 
tionship of V,H to V vc , simulation of the instrument response 
shows that angular features are resolvable for TJT\\ ratios 
greater than ~1.5). 

In Figure lb the response is for a bi-Maxwellian distribu- 
tion of 0 + ions with perpendicular thermal speed a factor of 5 
below the spacecraft speed ( V th -1 km s 1 and V xc = 5 km s 1 ). 
In the upper panel, even though V D = 0, the angular character- 
istics of the bi-Maxwellian are not resolvable within the de- 
tection thresholds of the instrument and the distribution ap- 
pears as a single peak at the spacecraft ram direction. The bot- 
tom panel assumes = 5 km s 1 parallel to B, the additional 




17,487 


GILES ET AL.: OUTFLOW IN THE AURORAL ZONE. POLAR CAP, AND CUSP 



drift speed of the distribution again apparent as a shift from 

In summary then, for plasma with thermal speeds much 
smaller than those of the spacecraft and drift speeds, the rela- 
drift motton dominates the response, leading to a 
rammed flux peak in the direction of relative wind as seen in 

uiven P h CeLra hH T 6 ' CaSe ^ angular info ™ation is 

t.ven by subtle features of the detailed shape of the rammed 

soeeds ,S h°[ ten ' ndetermina,e - For Plasma with thermal 
speeds much larger than those of the spacecraft and drift 

speeds, the- true plasma frame distribution function emerges 
clearly, although the relative drift motion, appearing as small 
asymmetries in the angular information, may be difficult to re- 
solve. The RIMS data included in this survey span this entire 
range ol extremes, although the specific set of plasma parame- 
ters contributing is mass dependent. 

Figure 2 illustrates this for H* and O* by plotting plasma pa- 
rameters covered by the data set as a function of ion Mach 
number. Here ,on Mach number is defined as the ratio of the 
d.str, button drift velocity (in the Earth frame) to the ion ther- 
mal speed (2 kT/M) ,n . Roughly speaking then, it is for sub- 
sonic (thermal speed greater than drift speed) ion Mach num- 
bers that the angular characteristics of the distribution will be 
resolvable^ Figure 2 shows that for drift velocities less than 

low Tc S | n ’ “ Utn LJ 0WS are subsonic for thermal components as 
s • e . Hence the angular characteristics of the H + 
samples in the data set should be easily resolvable except for 
the coldest plasma with strong (V D > 15 km s ') drift veloci- 
ties. O* outflows are subsonic at 15 km s' 1 for thermal compo- 
nents above 10 eV and at 8 km s- for thermal components 

i V | LV ™ erefore ’ for heav y ton outflows with tempera- 
tures below 5-10 eV, the angular information will frequently 
be indeterminate. M y 

Similarly, the drift motion effect is resolvable in the super- 
sonic Mach number range. Because for small thermal compo- 


nents ,t is the relative difference between drift and spacecraft 
velocity that is important. Figure 2 can be reinterpreted in 
terms ol a zero ion drift velocity. The shaded region then 
serves to delineate the range of spacecraft velocities with alti- 
tude (-7.6 km s ' at 1.5 R E geocentric, -4.2 km s ' at 3 R E , and 

7. km s at 4 5 and therefore the ion Mach range for a 
discernible relative drift motion effect to the spin angle mea- 
suremems (i.e., not dominated by the angular response) In 
the shaded region the relative difference between drift and 
spacecraft velocity is small, so that light ion (H<) outflow 
measurements with thermal energies below -0.4 eV will be 
dominated by the relative drift motion - a range for which de- 
tection is difficult due to positive spacecraft potentials For 
heavy ions (O*) the drift motion effect is significant for ther- 
mal energies up to 5 eV, and so will be a frequent spin angle 
feature in the heavy ion outflow samples. 

3.2. The Classification Scheme 

Our approach is to separate the spin angle data into cate- 
gories representing distinct physical behaviors, both for an- 
gular characteristics (the conventional pitch angle distribu- 
tion) and on the basis of apparent drift speed relative to the 
spacecraft. As a consequence to this, since the local ion es- 
cape speed is near to that of the spacecraft, the apparent drift 
speed allows us to categorize plasma behaviors loosely in re- 
1 at ion to the local escape speed. 

Figure 3 shows data measured by RIMS, and they are given 
as examples of the distribution types featured in this paper. 

I he data are plotted as accumulated counts per sample versus 
spin phase angle. In addition, the detector pitch angle is plot- 
ted as a function of spin angle. The counts per sample is pro- 
portional to the integral ion flux for the energy range 0 to 50 
eV with approximate conversion ratios of counts to integral 
lux of 31x10 coun,s c m 2 s ' ster' for H\ 1.9xl0 4 counts 
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DE-1 RIMS *1/287 20:37 - 20:38 DE-1 RIMS 82/071 19:34 - 19:35 



Spin Angle (Deg) Spin Angle (Deg) 

DE-1 RIMS 81/295 13:38 - 13:39 DE-1 RIMS 81/285 22:50 - 22:52 



Figure 3. Examples of DE l RIMS spin angle distributions typical of high-latitude outflow. Approximate 
conversion factors from count rate to integral ion flux are given in the text. Illustrated are: (a) unidirectional 
field-aligned beams, (b) uni-directional field-aligned conics, (c) asymmetric distributions characteristic of 
the upwelling ion signature, and (d) rammed distributions, which illustrate the effect of comparable space- 
craft and plasma outflow velocities as discussed in section 3. 


cnr 2 s l ster* 1 for He + , and 4.4x1 0 4 counts cm -2 s* 1 ster 1 for 0 + . 
These factors, which reflect the mass-dependent angular re- 
sponse, are based on preflight and in-flight calibration proce- 
dures; the ratios are time-dependent and are offered here only as 
a general approximation. 

The unidirectional field-aligned beam distributions (Figure 
3a) are characterized by a single flux peak certered at the field 
direction (here near 180° pitch angle). These distributions 
may be shifted from the field direction due to a combination of 
convection, drift, and spacecraft velocity effects. The impor- 
tant distinction is that these distributions are ordered by pitch 
angle (magnetic field direction) instead of by spin angle (the 
velocity vector of the spacecraft) and therefore correspond to 
field-aligned flows with velocity greater than that of the 
spacecraft. The rammed distribution (Figure 3d) is most com- 
monly identified in the data by a dome-shaped symmetric peak 
centered at or near the satellite ram direction. From the argu- 
ments above, the rammed distributions correspond to plasma 
flows that are slower than the spacecraft. The unidirectional 
field-aligned conic (Figure 3b) are those having flux peaks be- 
tween the directions perpendicular and parallel to the magnetic 
field direction. In this example, the spin angle distribution 
appears as a large single peak but with a distinct drop in count 
rate at the field direction (180° pitch angle). In other exam- 
ples, the conic can appear as two separate peaks. The asym- 
metric distributions (Figure 3c) appear as a single peak in 
count rate near the spacecraft ram direction but with a notable 
asymmetry toward the upward magnetic field direction. RIMS 
samples that have a skew coefficient (equation (5)) exceeding 
0.75 toward the upward magnetic field direction will be catego- 


rized asymmetric. In the example shown the coefficient of 
skew is 0.798 defined about the centroid position (14.6° in 
this example) for the portion of the distribution (61.7° in spin 
angle width) above the count level equal to 1/e 2 times the max- 
imum. It is this data signature that is associated with the up- 
welling ion events first identified by Lockwood et al. [1985a]. 
Distributions fall into the intermediate category when a single 
peak in the count rate falls between the field and ram direc- 
tions. For this reason, they cannot be readily assigned to ei- 
ther category and tend to be the samples in which the plasma 
and spacecraft velocities are more nearly equal. Within the 
RIMS data set, downflowing field-aligned beam distributions 
primarily consist of ions falling under the influence of grav- 
ity, the velocity of which is expected to be only a few kilome- 
ters per second. The unidirectional field-aligned downflows 
category is a subset of the rammed and intermediate categories 
and contains the ion downtlows for which a clear downward 
drift velocity is resolvable. 

The remaining distribution classifications covered in this 
paper, nothing observed and low counts, are to classify data 
signatures that do not fall into one of the categories defined 
above. The data intervals that have counts below the 2 counts- 
per-sample noise level are assigned to the 'nothing observed’’ 
category. The data intervals that have counts above the noise 
level, but without identifiable maxima are assigned to the 
"low-count" category. A later paper will examine distribution 
types primarily identified with midlatitudes (i.e., counter- 
streaming and trapped); all distribution types are included in 
the counting statistics so that the porportions reported are 
representative of the full population. 
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Figure 4 quantifies the inferred outflow drift velocities asso- 
ciated with three of the above categories by providing a rough 
estimate of the bulk drift velocity for each sample as a func- 
tion of spacecraft velocity (and consequently, as a function of 
altitude). The heavy black line marks the ion drift velocity 
necessary to escape Earth's gravitational potential. The inten- 
tion is to bracket the range of plasma drift speed that result in 
a given sample being assigned to one category or another. 
Other studies by DE investigators [Chandler et al . , 1991; 
Pollock el at., 1990] report detailed velocity statistics for po- 
lar wind and cleft ion fountain outflows. 

The calculation assumes the plasma drift velocity to be 
aligned with the magnetic field (i.e., no convection) and that 
only the spacecraft velocity component transverse to the field 
may bias the observed distribution toward the direction of 
spacecraft motion. We use the approximation that the mag- 
netic field lies in the spacecraft spin plane . This is generally 
valid within 12° for the present data set (and within 5° for 91% 
ot the samples). Provided the spin angle corresponding to the 
centroid of the observed plasma distribution (a) and the spin 
angle corresponding to the magnetic field direction (/J) do not 
coincide, the plasma drift velocity may be calculated from 

,z _ ^spacecraft s ' n ( a ) 

P ,asma sin(/J - a) (7) 

Figure 5 is a sketch of the vector measurements and angles 
relating to this calculation. For the rammed category all out- 
flow velocities, with few exceptions (all H + ), fall below the re- 
quired plasma escape velocity across the entire altitude range- 
consequently, the samples represented by this category, with- 
out additional energization, must either reenter the ionosphere 
or be trapped by magnetic mirror effects. The intermediate 
category is a mixture of escaping and nonescaping H + sam- 


ples, and contains mainly nonescaping 0 + . The beam cate- 
gory contains the majority of escape velocity samples for 
both ion species, although there is a significant portion of 
relatively 'fast," but nevertheless nonescaping, 0* at the 
highest altitude (lowest spacecraft velocity) bins. 

At high latitudes, then, the field-aligned outflows fall into 
three categories (rammed, intermediate, and field-aligned 
beams), which, for the energy range detectable by RIMS, cor- 
respond to ion flow speeds that are (1) less than, (2) compara- 
ble to, or (3) much faster than that of the spacecraft. This sep- 
aration by observed spin angle distribution serves to identify 



Figure 5. Sketch of the vector measurements and angles for 
the bulk drift velocity calculation shown in Figure 4. The 
value a is the centroid of the observed plasma distribution, 
and (i is the spin angle corresponding to the magnetic field di- 
rection. The calculation assumes the plasma drift velocity to 
be aligned with the magnetic field (i.e., no convection) and 
that only the spacecraft velocity component transverse to the 
field may bias the observed distribution toward the direction of 
spacecraft motion. 
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Geocentric Altitude (R £ ) Geocentric Altitude (R E ) 

Figure 6. (a) The percentage of samples, during which RIMS was programmed to accept a given ion, the 
instrument was able to detect that ion. Shown are seven Kp bins for H + , He + , 0 + , M/Z= 2, and 0 ++ . The occur- 
rences are calculated for invariant latitudes above 60° without regard to PAD category, altitude, or local time, 
(b) In six panels, occurrence with altitude for asymmetric, unidirectional field-aligned conics, unidirectional 
field-aligned beams, rammed, intermediate, and nothing observed; the calculations include only data above 
60° ILAT. The prominent feature to note is the separation of the various distribution types into distinct alti- 
tude ranges. 


which ions will be gravitationally bound and which are more 
energetic and thus able to escape to outer magnetospheric re- 
gions. 

4. Statistical Results 

4.1. Occurrences With Kp and Altitude 

The RIMS instrument, when programmed to detect a selected 
mass species, may or may not detect appreciable ion counts 
above the noise threshold within the 0-50 eV energy range. 
Figure 6a summarizes, as a percentage of the total number of 
samples, when RIMS measured ion counts above a 2 counts- 
per-sample threshold. Data are shown for seven magnetic ac- 
tivity level (Kp) bins for the ions, H\ He + , 0 + , M/Z- 2, and 
0 ++ . The occurrences are calculated for invariant latitudes 
above 60°, without regard to altitude, local time, or the cate- 
gories of pitch angle distributions defined above. The figure 
simply indicates in what percentage of samples, during which 


the instrument was programmed to accept a given ion, the in- 
strument was able to detect that ion. 

Counts of 0-50 eV H + are measured above the detection 
threshold in 75-78% of the samples, with no significant varia- 
tion with magnetic activity. In contrast, occurrences of the O* 
ions are strongly dependent on the activity level. At the low- 
est Kp bin, O* ions are measured in 33% of the samples; this 
percentage rises steadily with activity, exceeding the occur- 
rences of H + for the highest Kp bin. The He* ion occurrence 
probabilities lie between those of H + and O* and point to a 
weak dependence on activity. Displaying a pattern similar to 
that of O*, O** appears in about 19% of the samples for low Kp 
levels and increases to 41% during times of high Kp. M/Z= 2 
ions are seen in 20-29% of the samples, with no clear trend 
with activity. 

Figure 6b has six panels showing the percentage of occur- 
rences as a function of altitude for: (1) asymmetric distribu- 
tions, (2) unidirectional field-aligned conics, (3) unidirec- 
tional field-aligned beams, (4) rammed, (5) intermediate, and 
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MoTJt 8 °H SerVed; ' he Ca,Culati0 " s ' nc 'ude only data above 
60 ILAT and are sorted without regard to local time The 

Sme » y C T r a,l0n ‘ “ “™ d for e “ h ta " 

t ■ prominent feature to note in the altitude pro- 

de an e ly P es int ° distinct alti- 

tude ranges. The asymmetnc, intermediate, and rammed distri- 

bunons appear with greatest frequency at altitudes below 2.25 

m the h‘ sarnn. ‘ T * 0 " 31 conics - clearly identified 
175 and 7 ” Pr ° babi ' ity 31 ^idaltitudes between 
' 5 and ■; R f- Unidirectional beams tend to populate the 
h.gher end of the altitude range; they reach maximum fie 
quency above 2.5 * £ for H* and He* and 3.5 R E f or O T he 
noth.ng observed category, samples with counts below the 

0 P e akS Pr0bab " ity at bi 8 her altitudes for 

He , O , and O . Statistics for the M/Z=2 PAD types not eas 
■ y v-ewed in this format, display relatively higher frequenc es 
a. the lowest altitudes with the exception of intermediate 
distributions which occur with uniform frequency across the 
nt, re altitude range. It should be noted that in this format 
dts nbu, lons occurring in relatively small areas of local time 
or latitude will be averaged out in the sort over all local times 
and latitudes above 60‘; this is especially true for the M/Z-2 
'on in which the overwhelming number of nothing observed 
samples masks the contribution of auroral zone outflows. 
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4.2. Low-Altitude (R E < 3.0) Observations 
Plate 1 shows low-altitude (R E < 3.0) ILAT versus MLT oc- 

ZTJ/T abi JT S (2 5 ° by '' h0Ur bins) for (Plate 

la), field-aligned beams (Plate lb), asymmetric (Plate lc) and 
intermediate distributions (Plate Id). Results are shown for 
H He , and O and divided into Kp < 3 and Kp > 3. The linear 
scale a, the top of these plots and the ones that follow rep e 
sent the occurrence probability as a percentage, the number of 
samples luting the distribution criteria divided by the total 
number of samples for that ion. The color ,ed represents all 
probabilities greater than the number noted to the right- simi- 
larly^ the color blue represents all probabilities lower than the 
number noted to the left. Corresponding statistics for M/Z=2 

• w,,hout ><P division, are shown in Plate 2 Fewer 
M/Z=2 and O** samples at low altitudes preclude statistical cal- 
culations by activity level. Note that, for this altitude range 
the predawn local time quadrant between 60° and 75° ILAT is 
not well sampled ("white" indicates less than 15 samples were 
obtained in a given bin). Although the discussion applies 
only to the variability of distributions occurring above 60° 
LAT. we chose 50° ILAT for the minimum value in these spec- 
rograms to provide overlap and ease of comparison with mid- 
latitude statistics reported elsewhere. 

For R e < 3.0 there are relatively few occurrences of H* 
rammed (Plate la) at polar cap or auroral zone latiLII except 

»d Jl'n ,<M> !" ta«en W-W- 

1 l0n species display high occurrences for 

the rammed category below 60° to 65° ILAT. These are plas 
masphere isotropic distributions and characteristically Ld 
toward lower invariant latitudes with increased activity He* 
rammed PADs occur at greater than 50% probability from mid- 
night to dawn between 60° and 80° ILAT. The O* rammed oc- 
curs at >50% probability covering ILAT > 70° from pre 
m, mg t to about 9 MLT. In addition, the O* rammed occur' at 
>25% probability between 60° and 70° ILAT across the day- 
side. For all ions there is an absence of occurrences at auroral 
zone latitudes for afternoon local times; occurrences also de- 


crease with increasing magnetic activity. Statistics for 
rammed O ** fo" ow the O* pattern with the Lei of occurrence 
probability being similar. Af/Z = 2 follows the He* pattern 
with probabilities lower by -50%. In the case of M/Ll an d 
O and to an extent O*. lower probability at plasmasphere 
latitudes as compared to the H* and He* statistics will be re- 
flected as high probability for the same latitudes in the low 
counts category. This is a consequence of both the lower red 
live concentration of these ions within the plasma and instru- 
ment operational characteristics (less frequent sampling L 
suiting in erratic low-count distributions). 

H* low-altitude field-aligned beams (Plate lb) occur with 
grea er than 30% probability in an oval pattern that covers 
nearly a!l focal times and appears to be roughly coincident 
with the auroral oval. There is an expansion of this oval to 
lower latitudes for high Kp level and evidence of an increase in 
probab.ltty. The He* and O* probabilities also follow the au- 
roral oval pattern except for an absence of beams in the morn- 

ALn eC rh 0r ~ thC 3rea P ° PUlated by rammed distributions. 
Again, there is an expansion to lower latitudes for higher Kp 

and increased probability. O** statistics for the beam distribu 
t.ons are not conclusive, although it can be said that O** 

andTLaM Ur Wi,h a t b °L ,0% Pr ° bability at the same 'atiludes 

and 'oca! times as for O*. Since O* beam statistics increase 
niformly with activity, we might expect the averaging of O** 

pL a r2°TL al ^L t0 K l0wer the mean proportions - 

Plate 2. The M/Z= 2 beams, a category which also contains 

T“ ntS ° f ~ l keV protonSl follows th e auroral oval 
with -10% probability. 

For all ions the asymmetric PAD (Plate lc, and also referred 
to as upwelling ion events) has greatest occurrence on the 
prenoon dayside at invariant latitudes above 65° for the light 
tons and above 75° for the heavier ions, in addition, there are 
significant probabilities, increasing with Kp, in the afternoon 
aysi e and . for the light ions, in the region about 20-22 
MLT above 70° ILAT. The overall level of occurrence is 
higher during active conditions by as much as 20%; much of 

i, bei " 8 ‘ he ,ncreased afternoon dayside events 

Unlike H and He*, few O*. M/Z=2. or O** asymmetric events 
are observed for the polar cap. Lockwood e, al. [1985a] noted 
that when the dayside cleft upwelling feature is identified for 
O .the feature is present for all ions measured, so that the cleft 
region behavior exhibited by O* is representative of that for 
H . He M/Z=2, and O** as well. The reverse is not always 
true. When sorting the lighter ions, the programmed pattern 
recognition algorithm identifies additional periods of time as 
fitting the asymmetric criteria (skew coefficient >0.75) 
Therefore more samples and spatial coverage are tagged as 
asymmetric for H* and He* than for O*. The upwelling ion 
events associated with the dayside cleft region [Lockwood e, 

f; 198531 ™ clearly see " the O* statistics but are imbed- 
ded within a larger population of heat flux distributions for H* 
and He + . 

H* intermediate distributions (Plate Id) occur at 10-15% 
P 7 bab,l,ty across the dayside between 60° and 70° ILAT and at 
10-15% probability over the polar cap. He* and O* intermedi- 

65 L re 7 n° n'*i % Pf0bability ° Ver the P° lar ca P and between 
65 to 70 ILAT for all local times. Those over the polar cap 

increase in probability with activity. Statistics for O** fol- 
ows the O* pattern with a similar spatial pattern but somewhat 
lower probability. Few M/Z=2 intermediate events are ob- 
served above 55° ILAT. 

Plate 3 shows low-altitude (R f < 3.0) ILAT versus MLT oc- 
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currence probabilities for downflowing distributions (Plate 3a) 
and field-aligned conics (Plate 3b). Results are shown only for 
H + , He + , and 0 + and divided into Kp < 3 and Kp > 3. The 
downflowing distributions (Plate 3a) are difficult to distin- 
guish from the rammed category for the RIMS energy range 
because the downward drift speed of the ions tends to be rela- 
tively small. H + downflows are discernible at polar cap lati- 
tudes above 80° with an increase in occurrences with activity. 
There are also H + downflows at 70°-75° ILAT from 11 to 13 
MLT. CP downflows appear over the nightside polar cap above 
80° ILAT with no resolvable change with activity increase. 
Irregular regions of 0 + downflows also appear below 60° ILAT 
and are believed associated with convection near the plasma- 
pause. Occurrences of He + , A//Z=2, and CP+ downflows are less 
than 5% overall. The H + unidirectional conics (Plate 3b) are 
most frequently observed at invariant latitudes between 70° 
and 80° in the dayside local time sectors, and at invariant lati- 
tudes between 60° and 70° on the nightside. There is an in- 
crease (about 10%) in probability for observations of conics 
with an increase in magnetic activity. The unidirectional 
conic occurrences for the remaining ions are below 10% al- 
though, at least in the case of He + , the occurrences show hints 
of following the H + pattern. 

4.3. H ig he r- A I ti t ud e (3.0 < /?£ <.4.67) Observations 

Plate 4 shows high-altitude (3.0 < R E < 4.67) ILAT versus 
MLT occurrence probabilities for rammed (Plate 4a), field- 
aligned beams (Plate 4b), asymmetric (Plate 4c), and interme- 
diate distributions (Plate 4d). Results are for H+, He + , and 0 + 
and divided into Kp < 3 and Kp > 3. Corresponding statistics 
for M/Z— 2 and 0 ++ , without the Kp division, are shown in Plate 
2 . 

The high occurrences of heavy ion rammed distributions 
(Plate 4a) at low altitudes do not appear in the high-altitude 
statistics. In fact, there is almost a complete absence of 
rammed distributions at auroral zone and polar cap latitudes for 
all the ions measured with the exception of a small area of 
about 10-15% probability for high activity at 70°-80° ILAT 
postmidnight for H + and at 65°-75° ILAT premidnight for 0 + . 
Again, the probabilities below 60° to 65° ILAT are associated 
with plasmasphere isotropic distributions. 

H + high-altitude beams (Plate 4b) occur not only for auroral 
latitudes but at all latitudes greater than 70° across the full day- 
side sector. On the nightside, beams are observed with compa- 
rable probability only for auroral oval latitudes from midnight 
to dawn, the beams almost disappearing at from dusk to mid- 
night. He + beams occur with greatest probability (20-30%) at 
invariant latitudes between 60° and 80° between 6 to 19 MLT. 
There are somewhat lower, but still significant (15-18%), 
probabilities at 60° to 65° latitude from midnight to dawn. 0 + 
high-altitude beams occur with greatest probability during 
higher activity. These are at latitudes greater than 70°-75°, 
covering the entire polar cap between for all local times; and 
also from 60° to 65° between midnight and noon. For the 
high-altitude observations, in addition to an invariant latitude 
expansion with increasing magnetic activity, the beams are 
also observed at lower invariant latitudes for heavier masses. 
This expansion is between 5° and 10° for each of the mass 
species. For A//Z= 2, a category which also contains 
measurements of ~I keV protons, there is about 15% 
probability between 75° and 80° ILAT on the dayside and 
between 60° to 70° ILAT on the nightside. The probability for 





divided into Kp < 3 and Kp > 3. 
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unidirectional beams is less than 10% with the exception 
of a small area of about 15% probability at 70°-80° 1LAT and 
1 M3 MLT. 

The asymmetric distributions (Plate 4c) make up less than 
5% of the PAD samples at high altitude with the exception of a 
small area of about 10% probability for H* at 65°-75° 1LAT on 
the nightside for high magnetic activity. Similar to the ram 
and asymmetric results, the intermediate category (Plate 4d) is 
a much smaller proportion of the high-altitude sample popula- 
tion and, instead of appealing spread over the entire polar cap 
as for the low-altitude observations, they are limited to the 
nightside auroral zone between 20 MLT and 2-4 MLT. The dis- 
tribution of observations is similar for all ions with lighter 
ions appearing at slightly lower latitudes and movement to- 
ward higher latitudes with higher magnetic activity. 


4.4. Absence of Cold Plasma Observations 

Plate 5 shows low-altitude (R E < 3.0) ILAT versus MLT oc- 
currence probabilities for nothing observed (Plate 5a) and low 
count distributions (Plate 5b) and high-altitude (3.0 < R E < 
4.67) probabilities for nothing observed (Plate 5c) and low 
count distributions (Plate 5d). Results are shown for H + He + 
and 0 + and divided into Kp < 3 and Kp > 3. Statistics for 
M/Z-2 and 0 ++ , not sorted for activity, are shown in Plate 2. 

H low-altitude nothing observed events (Plate 5a) occur for 
latitudes above 70° ILAT across the nightside and above 80° 
ILAT across the dayside for low Kp and with lower frequency 
above 70° ILAT at all local times for high Kp. For high alti- 
tude the H* nothing observed events (Plate 5b) occur over the 
polar cap and, for low Kp , at auroral zone latitudes near mid- 
night. There is also significant probability below 60° ILAT at 
dusk, indicating a depletion of thermal ions in the plasma 
trough area there. He* low-altitude nothing observed events 
occur at >30% probability for all latitudes above 65° across 
the nightside and above 70° ILAT across the dayside. There is 
an expansion to slightly lower latitudes for higher Kp and 
slightly overall decreased probability. For high altitude the 
He + nothing observed events occur at >50% probability for all 
latitudes above 65°-70° ILAT at all local times. High-altitude 
°* not hing observed events occur at up to 75% probability 
over the entire polar cap and auroral zone area for low Kp\ for 
high Kp the CP nothing observed events drop dramatically in 
frequency and follow the auroral oval pattern. 0 ++ generally 
follows the pattern of CP, although without division with Kp it 
is not known whether a decrease in occurrence probability oc- 
curs with activity. A//Z= 2 follows that of He*. 

Low-count statistics are not significant at low altitude 
(Plate 5c), the observations that do occur appear at outer plas- 
masphere/plasma trough latitudes. For high altitude the low 
count category (Plate 5d) becomes important for the H* ion 
near midnight at auroral zone and polar cap latitudes. 

5. Discussion 

The low-energy (<50 eV) H\ He*, CP, M/Z= 2, and O** RIMS 
spin angle distributions for latitudes greater than 60° have 
been classified into the following pitch angle distribution cat- 
egories: unidirectional field-aligned beams, downflows, con- 
ics, asymmetric, rammed, intermediate, low counts, and noth- 
ing observed. 


5.1. Ionospheric Source Regions 

The statistical survey shows that low-energy outflow from 
the high-latitude ionosphere occurs for all the measured ions 
with a spatial and statistical distribution dependent on mag- 
netic activity. For convenience, we discuss the outflows in 

terms of three regions: the dayside cleft, auroral zone, and po- 
lar cap. K 

During quiet times, cleft region upwelling ion events (the 
asymmetric distributions in Plate I) are observed almost ex- 
clusively ,n the prenoon sector. The area of most probable 
observation becomes smaller and is bounded at higher invari- 
ant latitudes for the heavier mass species. These asymmetric, 
field-aligned spin distributions are identified as ions carrying 
an upward heat flux, with both parallel and perpendicular ve- 
ocity components, and so "upwell" in response to heating be- 
low the satellite [ Moore et at., 1985; Biddle el at 1 985* 
Lockwood et at.. 1985a], Moore e, at. [1985] demonstrated 
that the upwellings exhibit conic structure in the energetic tail 
of the distribution and suggested that the distinction between 
the unidirectional beam, conic, and upwelling categories may 
be a matter of the degree by which transverse acceleration af- 
fects the low-energy core of the distribution. The probability 
of observing upwelling events in the morning sector is rela- 
tively steady regardless of activity level which contrasts 
sharply with the large increase with activity seen with the au- 
roral 0 + outflowing beams (Plate lb). This is consistent with 
Lockwood et at. [1985a] in which these prenoon O* outflows 
from the polar cleft ionosphere were first identified as a major 
source ot heavy ions for the polar cap magnetosphere. There 
is, in addition, a distinct spatial change for cleft outflow with 
increased activity - at higher Kp levels upwelling events ap- 
pear on the afternoon dayside and the area of peak probability 
spreads to lower latitudes. In a study of the source characteris- 
tics of these events, Moore et at. [1986] found typical ener- 
gies to be of the order of 10 eV and found that the distinctive 
transverse ion heating exceeds gravitational escape require- 
ments. Evidence from this and the other studies indicates 
then, that the dayside prenoon cleft is a consistent source of 
low-energy ions of velocity sufficient to escape to outer mag- 
netosphere regions regardless of species (including the minor 
ions M/Z=2 and 0 ++ ) or activity level and that the dayside af- 
ternoon auroral zone becomes an additional source during in- 
creasing activity. 

While it is clear from the statistics that this asymmetric 
data signature, an indicator of transverse heating in the iono- 
spheric outflows, is specific to auroral latitudes, it cannot be 
unambiguously stated that they are unique to the dayside 
auroral zone. Figure 6b shows the peak observance of these 
events to be between 1.5 and 2.0 R E geocentric altitude. 
Unfortunately, due to a combination of instrument operation 
and orbital constraints, samples taken below 1.75 R E geocen- 
tric cover only the 5 to 12 MLT range. Samples between 1.75 
and 2.0 R E geocentric add auroral latitudes from 19 to 24 MLT. 

If, as has been suggested, the upwelling data signature evolves 
by magnetic folding to a field-aligned conic distribution, and 
if the asymmetric data signature does exist for the entire auro- 
ral zone, we would expect to find ample evidence of conics at 
nightside auroral latitudes, which is the case as can be seen in 
Plate 5b. 

Low-energy ion outflows associated with the auroral oval it- 
self are mainly of the unidirectional beam type (Plate lb). For 
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Plate 5. Low-altitude (R E < 3.0) ILAT versus MLT occurrence probabilities for (a) nothing observed and 
(c) low-count distributions and high-altitude (3.0 < R E < 4.67) probabilities for (b) nothing observed and (d) 
low-count distributions. Results arc shown for H\ He + , and Q + and divided into Kp < 3 and Kp > 3. 
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H , unidirectional beams are observed in the auroral zone for 
most all local times, increasing by about 10% with activity. 
The beam distributions are observed at progressively lower 
invariant latitudes for heavier mass species. These 
observations are in agreement with the earlier studies of 
Ghielmetti et al [1978] and Gorney et al [1981] in that the 
occurrence ol upward flowing beams appears to be associated 
with the statistical auroral oval. The occurrence frequencies 
peak more equatorward on the nightside than for the dayside, 
and peak at lower invariant latitudes during disturbed times. 
The observations are also consistent with Yau et al. [1984], 
who also tound that the oval of upflowing beams is energy de- 
pendent (the low-energy upflowing oval occurring poleward of 
the energetic ion upflowing oval). In addition, for heavier 
ions the beams display a distinct asymmetry in distribution 
speed between dawn and dusk local time sectors. For He + and 
O the unidirectional beams are observed primarily for the dusk 
sector and increase significantly in probability with increas- 
ing activity (as much as a 30% increase for O*). The "nothing 
observed" category complements this duskside behavior by 
having probability greater than 50% in this same region for 
low activity and decreasing sharply for high activity (Plate 
5a). In the dawn sector the heavier ion distributions, espe- 
cially O**. tend to be the lower velocity rammed type (Plate 
la). Since the unidirectional beam distributions generally 
have adequate escape velocity (Figure 4), the auroral zone then 
appears to be a steady source of escape velocity H + for all local 
times, a steady source of escape velocity He* ions on the 
duskside, and a source of escape velocity heavy ions on the 
duskside but primarily during increased activity. The fact that 
the He* and O* low-altitude dawn-sector rammed and 
intermediate distributions (Plates la and Id) fill in the missing 
beam distributions suggests that there is an amount of low- 
energy heavy ions unable to escape to high altitudes unless 
acted upon by energization processes which are correlated with 
increased activity. The M/Z=2 ion beam statistics contain 
measurements of -1 keV protons and so it is unclear whether 
escape velocity outflows other than the asymmetric type exist 
for this mass type. 

For latitudes poleward of the auroral zone the pitch angle 
distribution statistics vary with mass, activity level, and local 
time. The majority of the heavy ion observations in which 
significant counts are detected (>50 or 60%) are of the rammed 
and intermediate type (Plates la and Id), which still indicate 
field-aligned upflow, but with velocities comparable to the 
spacecraft and shown in section 3 not to be sufficient for es- 
cape to higher altitudes without additional energization. There 
is also a dawn-dusk asymmetry reflected in the ion distribution 
speed in which higher speed distributions appear with great- 
est frequency at dusk and lower speed distributions have great- 
est frequency at dawn. There appear to be small changes in the 
probability levels with magnetic activity, a slight increase in 
rammed distributions, and a small decrease in intermediate ob- 
servations; however, it is important to note that there is al- 
ways low-energy O* at low altitudes over the entire high- 
latitude region regardless of activity (Figure 6b). On the other 
hand, H + has a significant number of heat flux (asymmetric), 
unidirectional beams, and intermediate events for latitudes 
poleward of the auroral zone, the probabilities for these cate- 
gories increasing with activity. There are relatively few 
rammed distributions H*, although the "nothing observed" 
category for H* has >50% probability for low Kp. The He* ob- 
servations have characteristics of both H* and O*. the "noth- 


ing observed '/heat flux statistics of H* and the 
rammed/intermediate statistics of O*. Overall, the region 
poleward of the auroral zone is again a consistent source of 
low-energy ions, although we would expect only the lighter 
mass particles to have sufficient velocity to escape to higher 
altitudes. 

5,2. Evolution of Field - A ligned Distributions 
During Outflow 

By comparing distribution types at low altitudes with those 
at high altitudes we can infer something about how low-energy 
ion distributions evolve after ionosphere injection. The 
rammed distributions, which make up a significant portion of 
dawnside heavy ion observations at low altitude, essentially 
disappear at high altitude (Plate 4a). These either reenter the 
ionosphere or are trapped by magnetic mirror effects. If these 
low drift velocity rammed distributions instead experienced 
energization, such that the ramming effect no longer domi- 
nated, then we would expect a significant portion of "faster 
speed" distributions such as intermediate or field-aligned 
beams in the dawnside polar cap, which is not the case. For 
similar reasons the intermediate distributions seen at low alti- 
tudes are also believed to return to the ionosphere or remain 
trapped. At low altitude, only the asymmetric and field- 
aligned conic distributions, which reflect a significant degree 
of transverse acceleration, and the field-aligned beam distribu- 
tions appear to account for the statistical patterns observed at 
high altitudes. The observations support two mechanisms, 
which may well be related, for low-energy outflow into the 
outer magnetosphere (the cleft ion fountain, so named by 
Lockwood et al. [1985b], and its counterpart in other sectors 
of the auroral oval) what the present authors call the auroral 
ion fountain. 

The cleft ion fountain, as explained by Lockwood et al 
[1985b], consists of ionospheric plasma of 1-20 eV energy 
streaming upward into the magnetosphere where high-latitude 
convection electric fields can modify the ion path resulting in 
a poleward dispersion. This effect is a well-documented feature 
[cf. Waite et al. , 1986; Moore et al, 1985; Tsunoda et al . , 
1989] and is clearly evident by the beam distributions at high 
altitude over the polar cap (Plate 4b). It is not certain whether 
all the high-altitude O* polar cap observations originate from 
the cleft ion fountain mechanism; however, the very low 
probabilities of observing escape velocity O* outflow at low 
altitudes in the polar cap would indicate that this is so. The 
outflows in Plate 4b also demonstrate the mass spectrometer 
eltect associated with the cleft fountain whereby the lighter 
ions escape to the magnetosphere tail regions, but the heavier 
low-energy ions, such as O*, are more influenced by the con- 
vection electric field and are swept toward the nightside auroral 
region and plasma sheet. For weaker convection, heavier ions 
outflowing on the dayside would be swept antisunward over the 
polar cap and reenter the nightside polar ionosphere. This is 
supported by the Plate 3 O* downflows. 

The auroral ion fountain concept is reflected in the field- 
aligned beam and conic results. For the RIMS instrument the 
unidirectional field-aligned beams correspond to pitch angle 
distributions that have not undergone magnetic mirroring and 
can be represented by a drifting Maxwellian distribution with 
supersonic drift speed [Moore et al , 1985]. Traveling upward, 
these ions flow along high-latitude field lines that are open; or 
closed in a highly elongated fashion such that bouncing would 
occur over long time periods and be disrupted by nonadiabatic 
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processes. Low-altitude auroral latitudes are populated by 
beam distributions near dusk, "slower" intermediate or rammed 
distributions at dawn, and conics on the nightside (Plates 1 
and 3). High-altitude, late-aftemoon beams are almost absent; 
instead the region is represented by the intermediate and 
"nothing observed" categories (Plates 4 and 5). Tt is hypothe- 
sized that low-altitude field-aligned beams flow outward along 
auroral field lines, and for late afternoon, experience addi- 
tional acceleration, such that the ion energy distribution ex- 
ceeds the instrument detection range (>50-60 eV). Higher- 
energy beams, particularly those >1 keV, have been shown to 
have strong local time dependence, being concentrated near 
dusk at auroral latitudes [Gomey et ai , 1981; Yau et ai, 1984; 
Kondo et ai , 1990]. 

The sequence of events for both the cleft ion fountain and 
the auroral ion fountain is consistent with the theory that a 
continuum of behavior exists such that the upwelling distribu- 
tions evolve with altitude into conics and then high-altitude, 
field-aligned flows (cf. Figure 6). Conics correspond to up- 
ward field-aligned flows that have some additional perpendicu- 
lar velocity component. Conic distributions typical of the 
midaltitude nightside have been studied by Moore et ai [1985] 
and found to originate between 1000 and 2000 km altitude 
which is consistent with the altitude distribution shown in 
Figure 6. The study of Gomey et ai [1981], based on electro- 
static analyzer data at higher energies (0.09 to 3.9 keV), found 
that ion conics occur in the geocentric distance range of 1.4 up 
to the S3-3 apogee at 2.3 R Et or about 2500- to 8000-km alti- 
tude and that ion beams occur primarily above about 2.2 R E or 
7200-km altitude, with the frequency increasing with altitude. 
During quiet times the S3-3 conics were associated with the 
daytime polar cusp region, but were observed uniformly in lo- 
cal time during disturbed times. The older studies did not dis- 
tinguish between events peaked near 90° pitch angle (our 
asymmetric category) and those peaked between 120° and 170° 
pitch angle (our conic category). However, they are consis- 
tent with our observations in that our dayside asymmetric 
events exist for all activity levels, while those observations 
associated with our more narrowly defined conic category are 
found to increase with activity and are distributed in local 
time. The S3-3 results are also consistent in that quiet time 
ion beams have maximum frequency in the premidnight sector 
and are mainly confined to the dusk sector during disturbed 
times. 

5.3. Considerations on the Absence of Cold 
Plasma Observations 

The "low-count" and "nothing observed" categories are rep- 
resentative of two possible conditions: regions in which the 
ions have energy so low that they are repelled by positive 
spacecraft potentials or of regions where there are no ions 
within the detection range of the instrument, the low-count 
samples being those that occur near the limits of the instru- 
ment response (-0.1 ions cm 3 with no potential barrier). 
Analysis of RIMS data by Oisen [1989] established a density- 
potential relationship such that the satellite accumulates posi- 
tive charge on entering density regions below 1000 cm 3 , ris- 
ing slowly to about +1 V at 100 crrr\ and about +5 V at 10 
cm 3 . By preventing measurements of the coldest plasma 
components, this potential can mask out an isotropic 
background plasma or, when the potential is sufficiently high, 
can give the appearance that no cold plasma exists in a region. 


This effect is known to be a concern in the plasmapause and 
trough region, where densities are of the order of 7 to 70 cm’ 3 
and the spacecraft reaches potentials of +2 to +5 V [Olsen, 
1989]. At altitudes above the ionosphere in the polar regions, 
densities again can be sufficiently low so that observation of 
ions of classical polar wind speed is possible only during 
negative aperture bias instrument modes to offset the positive 
potential [Nagai et ai , 1984]. We cannot always assume, in 
low-density regions, that cold plasma does not exist just 
because it has not been measured. 

The RIMS H + nothing observed events at low altitude may 
be representative of an inability to observe the lowest-energy 
polar wind outflows in low-density regions. As detailed in the 
above discussion on the "auroral ion fountain," at high alti- 
tude near midnight, the H + nothing observed events are 
thought to be measurements in which the low-altitude outflows 
have been accelerated beyond the range of the instrument. The 
He 4 and O* events in the dusk auroral zone are probably a mix- 
ture of the two effects. 

6. Summary 

We have presented a statistical study of the occurrence fre- 
quency distribution of low-energy ion pitch angle distribu- 
tions in the near-Earth magnetosphere using measurements 
from the retarding ion mass spectrometer aboard Dynamics 
Explorer 1. A relationship is shown between pitch angle dis- 
tributions of different masses (H\ He + , 0\ M/Z- 2, and 0 ++ ) 
and velocity characteristics. The added spatial dimensions in 
this study relative to earlier work provide a more detailed pic- 
ture of ion transport events. Occurrence frequencies are deter- 
mined for all measured ions as a function of altitude, invariant 
latitude, magnetic local time, and magnetic activity ( Kp ). 
From this study, several characteristics are apparent. 

6.1. Overall Occurrence With Activity 

The probability of observing H + and Af/Z=2 distributions 
for 1LAT >60° was found to be independent of magnetic activ- 
ity. He + was found to be weakly dependent, while O* and 0 + * 
were strongly dependent on magnetic activity. The large in- 
crease in O* and 0 ++ observations with increasing activity was 
found to be composed primarily of escape velocity unidirec- 
tional field-aligned and upwelling distribution types. 

Upwelling (transversely heated, upflowing) ions occur pre- 
dominantly, but perhaps not exclusively, in the pre-noon sec- 
tor at 70°-80° ILAT. At higher Kp levels the events are ob- 
served with increased probability on the afternoon dayside and 
at slightly lower latitudes for all ion species. 

The occurrence of unidirectional beams, for all ion species 
including M/Z-2 , is consistent with the location and geomag- 
netic activity dependent expansion and contraction of the au- 
roral oval. The data suggest that beam distributions observed 
at higher altitudes are predominately associated with the up- 
welling distributions observed at lower altitudes and that the 
paths for ions injected at low altitude in the dayside auroral 
oval include some amount of convection antisunward over the 
polar cap, as suggested in the Lockwood et ai [1985a, b] cleft 
ion fountain sequence. 

Peak conic and upwelling ion observations occur predomi- 
nantly at geocentric distances of 1.75 and 2.0 R E , respec- 
tively, while field-aligned beams occur most frequently above 
about 2.0 to 2.25 R E , especially for the heavier ions. Since 
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angular distributions tend to become increasingly "beamlike" 
at higher altitudes because of magnetic focusing it is thought 
that the upwelling events may well evolve into the observed 
unidirectional conic and beam distributions as the ions travel 
outward along the magnetic field direction. 

6.2. Ionospheric Source Regions 

The dayside prenoon cleft is a consistent source of low- 
energy ions of escape velocity regardless of species or 
activity level and the dayside afternoon cleft, or auroral zone, 
region adds to this source for increased activity. 

The auroral oval as a whole appears to be a steady source of 
escape velocity H 4 ions at all local times, a steady source of 
escape velocity He 4 ions for the dusk sector, and a source of es- 
cape velocity heavy ions for dusk but primarily for increased 
activity. 

The region poleward of the auroral zone is a consistent 
source of low-energy ions, although only the lighter mass par- 
ticles appear to have sufficient velocity to escape to higher al- 
titudes. 

6.3. Evolution of Field-Aligned Distributions 
During Outflow 

The observations support two concepts for outflow: (1) The 
cleft ion fountain consists of ionospheric plasma of 1-20 eV 
energy streaming upward into the magnetosphere where high- 
latitude convection electric fields modify the ion paths result- 
ing in a poleward dispersion of the ion trajectories. (2) The 
auroral ion fountain consists of field-aligned beams which 
stream outward along auroral latitude field lines, and in addi- 
tion, for late afternoon local times, experience additional ac- 
celeration, such that ion energy distributions tend to exceed 
the detection range of the instrument (>50-60 eV). 

6.4. Absence of Cold Plasma Observations 

The low-count" and "nothing observed" trends with mag- 
netic activity appear to complement that of the unidirectional 
field-aligned beam distributions and are representative of two 
possible conditions: regions in which the ions have energies 
relative to the spacecraft so low that they are repelled by posi- 
tive spacecraft potentials or regions with no ions having en- 
ergies within the detection range of the instrument. We can- 
not always assume, in low-density regions, that cold plasma 
does not exist just because it has not been measured. 
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The dynamics of earth’s plasmasphere is affected by gradients in the plasma pressure, the convec- 
tion electric field, and the corotation electric field. The Retarding Ion Mass Spectrometer (RIMS) on 
Dynamics Explorer 1 measures the pitch angles of this plasma from which the total plasma velocity 
vector can be determined. Using assumptions concerning the relative importance of convection and 
corotation, the velocity of ions along the magnetic field line can be determined. These results repre- 
sent the first direct observations of ion velocities along magnetic field lines at high altitudes in the 
plasmasphere. Data from October 1981 through January 1982 have revealed several interesting aspects 
of the field-aligned flow of H + and He + . Downward directed velocities observed in the northern 
hemisphere within 30° of the equatorial plane are interpreted as interhemispheric flows. Typical 
velocities are of the order of a few hundred meters per second. Higher-speed (up to 1 km s"l) flows 
observed in the outer plasmasphere appear to be associated with refilling. Counterstreaming of H + and 
He + has been observed in the autumn and winter hemisphere with He + flowing upward and H + flow- 
ing downward. Differences between these results and theoretical work can probably be attributed to 
differences in the season and the density level within the flux tube. 


1. INTRODUCTION 

The plasma in earth’s plasmasphere is generally Max- 
wellian with energies of less than a few electron volts. In 
order to adequately describe this plasma it is necessary to 
know several basic quantities, namely, (1) composition, 
(2) number density, (3) particle flux, (4) temperature, and 
(5) heat flux. These quantities have been measured in the 
ionosphere using incoherent scatter radar and low-altitude 
satellites. The composition, density, and temperature have 
also been measured in the plasmasphere by numerous 
satellites [e.g., Serbu and Maier, 1970; Decreau et ai, 1982; 
Chappell et al , 1971; Chappell, 1982; Comfort et ai, 
1985], but the ion velocity has not been measured exten- 
sively at altitudes above the ionosphere. While ionospheric 
observations have made important contributions [e.g., 
Behnke and Harper , 1973; Behnke and Kohl, 1974; Vickrey 
et al, 1976], they have been limited to altitudes typically 
less than 1000 km which is normally below the 0 + -H + 
transition height. Chemical processes are still important at 
and above the 0+-H+ transition height which can sig- 
nificantly affect the flow velocity for H + and 0 + and make 
it difficult to infer the velocity or flux at altitudes above 
the ionosphere. Using a one-dimensional model of the 
plasmasphere, Richards et al [1983] have shown that ion 
flow reversals can occur at altitudes just above the iono- 
sphere which can lead, for example, to incorrect estimates 
of refilling rates based on ionospheric observations. There- 
fore, to accurately determine such parameters as the rate 
of tube refilling or of interhemispheric transport it is 
necessary to observe the ion flux at high altitudes and in 
the equatorial region. 

Models of the plasmasphere have predicted many mter- 
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esting and significant features of the magnetic field-aligned 
plasma flow in the plasmasphere. The flow of plasma into 
and out of the equatorial reservoir has been modeled 
extensively [e.g., Banks et al, 1971; Moffett and Murphy, 
1973; Marubashi and Grebowsky , 1976; Murphy et al, 
1976]. In general, these models indicate that when the flux 
tube is significantly depleted of plasma, the ions (primarily 
H+) flow upward from the ionosphere into the plasma- 
sphere at all local times. As the flux tube fills there is a 
transition to downward flow at night. The modeled H + flux 
is typically found to be about 10^ cm’2 $'l at 1000 km and 
as large as a few times 10? cm _ 2 at high altitudes 
(>4500 km), corresponding to velocities of a few hundred 
meters per second. Murphy et al [1976] suggested that 
these flux values vary little as the tube fills. Banks et al 
[1971] and, more recently, Khazanov et al [1984] have 
modeled the time dependent nature of refilling and found 
that after about 2 days the H + velocity along most of an 
/,= 4 flux tube was subsonic. Above about 5000 km the 
Khazanov et al. model predicted H + velocities in the range 
from 0.1 to 1.0 km s’l and a flux of about 5x10? cm~2 
s‘l at 5000 km, consistent with the results of Murphy et a l 
[1976]. 

Models have also shown that the ion velocities are not 
generally symmetric about the equatorial plane and sig- 
nificant interhemispheric flows can occur [e.g., Richards 
et al, 1983], Such uneven flow patterns are the result of 
asymmetric illumination of the conjugate regions of the 
ionosphere, the seasonal variations in neutral composition 
and concentrations (e.g., the winter anomaly and the 
winter helium bulge), and the neutral wind pattern. The 
interhemispheric flow of H+ and 0 + has been modeled 
with typical values for the H + flux being of the order of 
106 to 107 cm'2 s'l, corresponding to velocities of a few 
tens of meters per second [e.g., Bailey et al., 1978; 
Chandler et ai, 1983]. The flux of He+ obtained by 
Richards et ai [ 1983] was more than an order of 
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Fig. 1. Illustration of the components of the total ion velocity as seen by RIMS. 


magnitude lower than the H+ flux, while the He+ velocity 
was, at times, higher. In addition, they showed that at 
times the flows of H+ and He + were in opposite directions 
(counterstreaming). 

Vickrey et al [1976] observed counterstreaming of 
0 + and H + at ionospheric altitudes above Arecibo. Bailey 
et al. [1977] modeled the results of Vickrey et al. and 
found counterstreaming primarily at night. Later, Bailey 
[1980] showed that counterstreaming should exist during 
most of the day over a range of altitudes from just above 
the F 2 peak up to the 0+-H+ transition height. These 
results were attributed to the flow of 0 + upward into a 
chemical sink in the transition region (due to resonant 
charge exchange with H) and the subsequent flow of H + 
downward from this chemical source. Young et al. [1979] 
showed that 0 + -H + counterstreaming should occur, in the 
steady state situation, during conditions of field tube 
symmetry and unequal ion and neutral temperatures. 

The results of Richards et al. [1983] suggest that 
H + -He + counterstreaming should exist over large portions 
of midlatitude flux tubes under a variety of conditions. 
In their simulations of refilling, H + flowed upward from 
both hemispheres into the plasmasphere, while He + flowed 
predominantly from the winter hemisphere as a result of 
the winter helium bulge. This produced a region of H + -He + 
counterstreaming in the summer hemisphere. The velocity 
of both ions was of the order of a few tens of meters per 
second in the region of counterstreaming. 

While the velocity of thermal plasma along magnetic 
field lines has been modeled extensively, it has never been 
measured in detail at altitudes above the topside iono- 
sphere. Using data from the Retarding Ion Mass Spec- 
trometer (RIMS) onboard the Dynamics Explorer 1 satel- 
lite, it is now possible to correct this deficiency. This report 
outlines a technique for determining the field-aligned 
velocity of 0+, H+, He + , He++, 0++ N+ in the plasma- 
sphere. In addition, results are presented from observations 
of H + and He+ in the high-altitude plasmasphere. These 
data are taken from the first 4 months of RIMS observa- 
tions (October 1981 to January 1982) and cover magnetic 
local times from 0600 to 1000 hours and from 1500 to 
1900 hours in the period between autumnal equinox and 
winter solstice. 


2. TECHNIQUE 

The technique used is based on the geometry of the DE 
1 orbit (specifically the relationship of the satellite velocity 
vector to the magnetic field vector). The DE 1 orbit plane 
is very nearly (<10°) parallel to planes of constant mag- 
netic longitude. Thus the RIMS radial head samples almost 
all pitch angles in each 6-s spin (see Chappell et al. [ 1981 ] 
for a complete description of the RIMS instrument). In 
order to measure the naturally occurring field-aligned flow 
velocity, it is necessary to first determine the direction of 
the total velocity vector, £tot> the plasma in the spin 
plane (see Figures 1 and 2). This is done by finding the 
average pitch angle of the plasma or, equivalently, by 
finding the average spin angle of the plasma distribution 
measured by RIMS. Subtracting the satellite ram velocity 
from this vector will leave, in general, a combination of 
flow vectors resulting from three separate mechanisms: 
corotation, convection, and magnetic field-aligned flow. 
Alternatively, by determining the magnitude and direc- 
tion of ^toti both the field-aligned velocity and the con- 
vection velocity could be determined. However, in practice 
finding the magnitude of ^tot is a time-consuming process 
requiring detailed fitting of theoretical spin curves to the 
data. In order to avoid this large increase in processing time 
a convection velocity is assumed and ^tot ls determined 
from the data. The method for determining ^tot anc * f° r 
separating the field-aligned velocity from the other com- 
ponents is outlined in the appendix. 

If £xo\ is taken as unknown, some assumptions must be 
made concerning the corotation and convection velocities 
in the absence of electric field data. The corotation velocity 
of the plasma inside the plasmasphere varies as a function 
of latitude and radial distance from earth’s center and can 
be as large as 2 km s _1 within the radial distance covered 
by the DE 1 orbit. However, the geometry of the DE 1 
orbit is such that in the plasmasphere the satellite travels 
along a constant magnetic longitude, to within a few 
degrees. Therefore plasma flow due to corotation will have 
no significant effect on the measured flow vector in the 
spacecraft x-y plane. 

Corotation will, however, have an effect on the 
estimated flux values. As DE 1 moves to larger radial dis- 
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Fig. 2. The Dynamics Explorer 1 spacecraft coordinate system. 


tances, the corotation velocity becomes comparable to the 
spacecraft velocity, and, as a consequence, the ions come in 
at an angle to the instrument central plane (the satellite 
spin plane). As a result of the decrease in response of the 
instrument as the angle increases away from the central 
plane, the measured count cate will decrease, and hence the 
estimated flux will decrease. For the types of orbits repre- 
sented in this paper, the decrease in instrument sensitivity 
will be less than 50%, or, in other words, the flux will be 
underestimated by no more than a factor of 2. This differ- 
ence is not important to the results of this study. 

Unlike corotation, convection velocities resulting from 
east-west electric fields are primarily in the spin plane of 
the radial detector. As such, their contribution to ^tot 
cannot generally be separated from that of the field-aligned 
velocities by the technique outlined in the preceding equa- 
tions. However, convection velocities in the plasmasphere, 
particularly the inner plasmasphere, are typically small, 
being of the order of a few tens of meters per second at F 
region heights [Richmond et al, 1980], Maynard et al. 
[1983] have shown that the plasma is primarily flowing in 
the direction of corotation in the plasmasphere for L< 4. 
Furthermore, if abnormally large convection velocities are 
anticipated, by restricting the analysis to portions of the 
orbit where the angle between the spacecraft velocity and 
the magnetic field approaches 90°, the influence of any 
radial convection velocity can be minimized. This results 
from the fact that near 90°, convection will be in the 
direction of the satellite velocity and as such will represent 
a ‘‘correction” to the ram velocity. In the plasmasphere 
the ram velocity is between 3 km s"l and 8 km s'*, while 
typical convection velocities are less than 1 km s - l [e.g., 
Maynard et al, 1983]. Therefore, the “correction” will 
be small. Thus for most work in the plasmasphere, convec- 
tion can be ignored, and it has been ignored for all results 
presented in this paper. 


3. Error Analysis 

The uncertainty in the determination of ^tot ant ^^B 
is a function of the uncertainty in all the parameters that 
enter into equation (14) (see the appendix). In order to 
estimate the typical errors involved in the analysis, use 
was made of a computer program which simulates the 
RIMS observations in the plasmasphere and generates data 
files identical to those created in the standard RIMS data 
analysis procedures. This allowed for the estimation of the 
errors associated with the numerical integration technique 
and provided a means of verifying that the analysis code 
was performing properly. The program simulates the output 
from a retarding potential analyzer for a given set of 
spacecraft-specific parameters (e.g., velocity, potential) 
and plasma-specific parameters (e.g., temperature, density, 
flow velocity) under the assumption that the plasma is 
Maxwellian. It was modified to allow for a total ion velo- 
city vector with a direction different from the ram direc- 
tion. 

To estimate the errors from the numerical integration, 
synthesized data were analyzed for the case where all 
spacecraft and plasma parameters were held fixed. Using 
1 hour of simulated data and 1-min averages of the “data,” 
the results for three different total velocity vectors show a 
typical error in determining the spin phase angle of the 
total velocity vector of about 0.15°. In addition to these 
ideal data, which produce smooth spin curves, synthesized 
data which included fluctuations in the flux as a function 
of spin angle were analyzed. The technique was able to 
determine the velocity direction to an accuracy of between 
0.5° and 1.0°. This translates to an estimated error in the 
velocities of about 100 m s'*. 

The indication from these results is that the integration 
error is constant with respect to the magnitude of the 
velocity. This is as expected since the uncertainty is a 
function of the second derivative of the spin curve with 
respect to spin angle, for a two-point finite difference inte- 
gration technique, and also a function of the integration 
step size in spin angle (taken to be 2°). The uneven spin 
angle coverage of the instrument during the accumulation 
period and the subsequent averaging into 2° bins also con- 
tributes to the integration error. Because the factors contri- 
buting to these errors are independent of the plasma 
parameters, the 0.15° uncertainty should be the minimum 
error. In addition, as expected for the case where the only 
error is in determining 0 O , the fractional uncertainty in 
the field-aligned velocity derived from these synthesized 
data is approximately equal to the fractional uncertainty 
in 0 O . For cases where the angle between the magnetic 
field vector and the spacecraft velocity was allowed to vary 
this equality in the uncertainties held for all angles. This 
relationship is understandable since the velocity is basically 
a function of the tangent of the angle between the ram and 
the total velocity vector which can be assumed to be linear 
for these small angles (<5°). A full error propagation 
analysis shows that the error in 0 O is by far the largest 
contributor to errors in the determination of the velocity. 
Thus an estimated error of 100 m s"l is assumed for all 
results. The magnetic field data provided in the spacecraft 
data are accurate to 1.5 nT [ Farthing et al, 1981] or less 
than 1% for typical plasmaspheric magnetic field magni- 
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Fig. 3. The orbital configuration for October 8, 1981, November 29, 1981, and January 25, 1982. 


tudes. The attitude of the spacecraft is known to within 
±0.25° [Hoffman et a!,, 1981]. 

There are additional uncertainties in the derived values 
,£tot and which arise from the correctness of the 
assumptions concerning the corotation velocity and the 
convection velocity. These uncertainties are more difficult 
to estimate with RIMS data alone. However, as more DC 
electric field data from DE 1 become available, these effects 
can be accurately included in the analysis. By restricting 
the analysis to cases where the angle between the orbit 
plane and the longitudinal plane of the magnetic field is 
small, the effects of corotation can be minimized. In fact, 
calculating the corotation Velocity, assuming strict corota- 
tion, and projecting this into the spacecraft x-y plane for 
the case where the angle between the orbit plane and the 
longitudinal plane of B is between 12° and 6° yields 
components of less than 10 m s _ l in the field-aligned 
direction, which is much less than the estimated overall 
error of 100 m s'l. 

The effects of convection velocities on the analysis are 
harder to estimate. However, by restricting the analysis to 
cases where the spacecraft velocity vector is sufficiently 
far away from the field-aligned direction these effects can 
also be minimized to the point where they can be ignored 
at L values in the plasmasphere ( L<6 ). In addition, most 
of the observations shown in this paper were made near 
dawn and dusk where the convection velocity is very nearly 
in the direction of corotation [cf. Maynard et al., 1983]. 

4. Observations 

The DE 1 satellite was launched into an elliptical, polar 

orbit with the orbit plane in the 1000-2200 LT meridian. 
This plane precessed toward earlier local times. The 
observations shown in this section are from the first 4 
months of the mission and cover local times from 0600 to 

1000 hours and 1500 to 1900 hours. Figure 3 shows the 
orbital configuration for October 8, 1981 (dayside), 
November 29, 1981 (dayside and nightside), and January 
15, 1982 (afternoon) which illustrates the extremes of the 
orbits presented here. The dayside passes through the 
plasmasphere range from 1=2.0 to L=5. 5 and from 0° to 
37 north magnetic latitudes, while the nightside passes 


range from L=2A to 5.5 and 28° to 56° north magnetic 
latitude. 

In the following results the upper limit in L represents 
the limit of the observed rammed-type plasma distribution, 
marked in the figures by a cross. When a sharp density 
gradient is present, its position is also noted. When no 
sharp density gradient is present, the plasmasphere is 
assumed to extend to the limit of the observed rammed 
plasma. The lower L limit is set where the detector satura- 
tion effects reach about 25%. This is done to avoid any 
possible bias in the spin distribution due to detector satura- 
tion. An example of the observed flow velocities is shown 
in Figure 4a. This is a pass through the dayside plasma- 
sphere from />=2.0 to />= 2.8 (i.e. f the location of the 
plasmapause as defined by the sharp density gradient with 
£). Positive velocities represent south-to-north flow parallel 
to B in all results presented in this paper. The velocities 
range from 600 m s~l to <100 m s~l and show a general 
decrease with increasing L and latitude. The H+ velocity is 
higher than the He + velocity, which is a common feature in 
most of the cases that have been investigated. The H + 
flux is about 10^ cm*2 S -1 in the inner plasmasphere, falling 
to about 10 4 * * 7 cm-2 S -1 in the outer plasmasphere. The He+ 
flux is about an order of magnitude less in both regions. 
In addition to the pass in Figure 4a , the results from two 
subsequent passes through the dayside plasmasphere are 
presented in Figures 4b and 4c. These three passes are 
representative of the midmorning sector, after the plasma- 
sphere has had several hours to respond to sunrise. Note 
that in Figure 4b , velocities have been calculated, and 
shown, for the rammed plasma observed outside of the 
location of the plasmapause density gradient (marked by 
the arrow). Rammed plasma was also observed outside of 
the location of the density gradient in Figures 4a and 4c, 
but the count rates were too low to provide reliable 
velocity determinations. 

The next group of passes, in Figure 5, are representative 
of the sunrise period and are from late November 1981. 
The DE orbit had precessed, by this time, so that it passed 
through the dayside plasmasphere at local times near 0600 
hours. Only the southern hemisphere was totally sunlit at 
the time of these observations with a shadow height in the 
northern hemisphere, for these flux tubes, ranging from 
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Fig. 4. Field-aligned velocities for H + (open squaies) and He + (solid circles) as a function of L for October 8, 1981. 
Positive velocities are from south to north. The arrows indicate the location of the plasmapause density gradient, and 
the crosses mark the outer boundary of the observed rammed distribution. 


100 km on November 25 to 350 km on November 29. 
Influenced by the asymmetrical illumination of the iono- 
sphere, an upward flow from the southern hemisphere over 
into the northern hemisphere would be expected. Indeed, 
a south-to-north flow was observed between 4° and 35° 
magnetic latitude, which is interpreted as interhemispheric 
transport. In addition, the velocities were generally higher 
than those observed at 0900-1000 LT (Figure 4). These 
data point out some of the different flow characteristics 
observed. Typically, at low L values the H -1- velocity is 
observed to be just slightly higher than the He+ velocity. 
However, in the outer plasmasphere the situation varies, 
depending on the condition of the plasmasphere. For 
example, in Figure 5 a this near equality in the two 
velocities holds out to the plasmapause. In contrast, in 
Figure 5c the H+ velocity is generally higher than the He + 
velocity between L - 4 and L =5, while in Figure 5d , He + is 
flowing faster than H + in approximately this same L value 
range. 

The three passes in Figure 6 all have the same UT-LT 
separation, and, if strict corotation is assumed to hold, they 
should represent three observations of the same region of 
the plasmasphere, each separated from the other by 2 days. 
These results exhibit the time variation of the ion velocity 
as a function of L . For example. Figure 6 a shows a case 
where the velocities are low and fairly constant as a 


function of L. In Figure 6b the situation is substantially 
different (probably as a result of a large magnetic storm 
which occurred between these two observations) with a 
flow reversal observed near the plasmapause. In contrast, 
Figure 6c shows an flow, with speeds up to 1 km s"l, 
passing through an almost stationary (or perhaps a counter- 
streaming) He+ population. 

The previous examples are taken from the dayside 
plasmasphere. During the time of these observations, DE 1 
also sampled the nightside plasmasphere at local times near 
1800 LT. Analysis of the data from the nightside revealed 
some interesting cases of H + -He + counterstreaming. Some 
examples are shown in Figure 7. These observations were 
made at mid-latitudes in the northern (autumn/winter) 
hemisphere, The observations in Figures la and lb were 
made during the postsunset period, while Figure 7c is from 
the midafternoon period. In all cases, He + is flowing 
upward out of the northern hemisphere, while H + is 
flowing downward in the northern hemisphere. 

5. DISCUSSION 

The results presented in section 4 were taken from a 
data base of approximately 50 such passes through the 
plasmasphere and were chosen specifically to illustrate 
certain phenomena, some of which have been theoretically 
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Fig. 5. Velocity of H + and He + as a function of L for (a) November 25, 1981; (b) November 26, 1981; (c) November 
28, 1981; and (d) November 29, 1981. 


predicted. These results do not constitute a large enough 
sample from which to draw general, statistically significant 
conclusions. However, interpreting these observations in the 
context of such well-studied plasmaspheric parameters as 
the density content of the flux tubes and the plasmapause 
position should provide a familiar basis from which to point 
out the importance of ion velocity observations in under- 
standing the dynamic nature of the plasmasphere. Addi- 
tionally, the variations of the ion velocities within the 
plasmasphere must be understood in terms of the variations 
in the plasmasphere itself resulting from such influences as 
the diurnal cycle, the time history of magnetic activity, 
and seasonal effects. 

Convection and the Diurnal Cycle 

During quiet magnetic periods the plasmasphere coro- 
tates with earth, while, under active magnetic conditions, 
convection compresses the plasma within nightside flux 
tubes and pushes it downward into the underlying iono- 
sphere and plasma on the dayside is convected away on 
paths which intersect the magnetopause, resulting in the 
loss of plasma [e.g., Park , 1970], The characteristics of 
the plasmasphere and the position of the plasmapause 
are, as suggested by Carpenter et al [ 1969] and confirmed 
by Chappell et al [1971 ], related to the level of magnetic 
activity present during the passage of a given flux tube 


through the night sector (2200-0600 LT). Variations in 
the strength of convection on time scales shorter than 24 
hours can produce local time asymmetries in the shape 
and content of the plasmasphere. 

An example of local time variations in the shape of the 
plasmasphere can be seen in the data in Figure 4. While 
these three passes cover the same local time, latitude, and 
altitude regions, they occur at different magnetic longitudes. 
It is immediately obvious that distinct differences existed in 
the morphology and dynamics of the plasmasphere as a 
function of longitude on this day. First, the plasmapause 
location is different for Figure 4 c as compared to Figures 
4a and 4b. In all cases the classical plasmapause, defined by 
a sharp density gradient with L , was evident by a drop in 
the count rate (not shown) by more than 2 orders of magni- 
tude in less than 0.25Z. However, in all cases there existed 
rammed plasma in the region just outside of this density 
gradient, and this is interpreted as the newly formed, outer 
plasmasphere in which refilling was taking place, similar 
to the results presented by Horwitz et al [ 1 984] . However, 
only in Figure 4b were the count rates high enough to give 
reliable velocities. 

As further evidence of the local time differences, notice 
that while the velocities are comparable in Figures 4a and 
4b (excluding the region of refilling in Figure 4 b), they are 
substantially greater in Figure 4c, particularly in the region 
outside L- 2.5. These differences cannot be attributed to 
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Fig. 6. Ion velocities for H + and He + for (a) October 13, 1981; ( b ) October 15, 1981; and (c) October 17, 1981. 


the current level of magnetic activity as measured by Kp 
(tor these three passes Kp was 4-, 4+, and 4-, respectively) 
but can be related to the level of magnetic activity during 
the passage of these flux tubes through the nighttime sector 
(2200-0600 LT). The Kp history over the past 4 days 
(see Figure 8; note that the 2200-0600 LT period is 
marked by a black bar) suggests that the three regions 
recently experienced significantly different modifications 
due to convection associated with high magnetic activity. 
On October 5, more than 2 days prior to the observations 
in Figure 4a, the daily sum of Kp was 13, while on October 
6 the daily sum was 12-. However, as the region in Figure 
4 Q rotated through the afternoon sector on the day prior 
to the observations (October 7), Kp had increased to 4, 
resulting in some depletion of the outer plasmasphere in 
‘ius local time region as plasma was convected outward. 
As this region rotated through the nightside, the magnetic 
activity increased further (to A'p=6+), compressing the 
plasmasphere inward on the nightside resulting in move- 
ment of the plasmapause to A = 2.8. The plasma region in 
Figure 4b experienced similar activity levels as a fucntion 
of local time and had approximately the same plasmapause 
location as that in Figure 4a. In contrast, the region in 
Pig ure 4c experienced similar absolute levels of activity (as 


measured by Kp) but at difterent local times than those in 
Figures 4a and 4b. The region in Figure 4c, on the day 
before these observations, was influenced by increased out- 
ward convection on the dayside, probably resulting in some 
plasma loss. Unlike those in Figures 4a and 4b, as this 
region passed through the nightside sector just prior to the 
observations, Kp was decreasing, resulting in less compres- 
sion and a more extended plasmasphere. In addition, the 
region was near local dawn on October 7, when Kp reached 
6+, and, as shown by Chappell et al. [1971], the plasm- 
pause position is little affected by high magnetic activity at 
these local times. On the other hand, it should have exper- 
ienced more outward convection and plasma loss on the 
dayside than the two other regions. 

The velocity versus L profile in Figure 4c is significantly 
different from that in Figures 4 a and 4b. The velocities in 
Figure 4c are higher throughout the pass than those in 
Figure 4 a or 4b and the He+ velocity increases with increas- 
ing L between L- 2.5 and 3.1. While the higher velocities 
might indicate refilling in this region, the fact that the H + 
and He+ velocities are essentially the same suggests some 
process not strongly related to the ionospheric sources 
(which are quite different for the two ions). These types of 
flows may be indicative of redistribution of plasma within 
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Fig. 7. Ion velocities for H + and He + for (a) November 29, 1981 (1900 UT); ( b ) November 29. 1981 (1900 UT);and 
(c) January 25, 1982 (1500 UT). 

the tube at altitudes above the ionosphere. For example, if the influence of the diurnal cycle should have become more 
the plasma within a flux tube is convected outward so that important. That is, as plasma flows upward from the iono- 
it resides in a tube of greater volume, density irregularities sphere following sunrise, to replenish the plasma lost during 
may occur along the flux tube [e.g., Richards and Cole , the night, higher velocities would be expected and they 
1979]. These irregularities would then dissipate, through should be directed from the summer hemisphere which is 
flows, as the plasma assumes a new equilibrium distribu- illuminated first. In Figure 5 the velocities for four passes 
tion. through the plasmasphere at local times near 0600 are 

As the DE 1 orbit precessed to local times near sunrise, shown. In contrast to the results of Figure 4 (for 1000 LT), 
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Fig. 8. The Kp history for the regions of the plasmasphere observed in Figure 4. The time lines at the top of the figure 
give the correspondence between local time and universal time (bottom axis) for the particular segment of the plasma- 
sphere sampled by each pass assuming strict corotation. The times are calculated backward in time using the LT-UT 
difference at the time of observation. 


the velocities are generally higher throughout the region 
observed. In addition, He"*" was observed to be flowing at 
approximately the same speed as H+ at the lower! values, 
in contrast to the 1000 LT observations in which the He + 
velocity was measurably lower than the H + velocity. Both 
of the results suggest a more dynamic plasmasphere at 
sunrise, as expected. 

Because the orbital period of DE 1 is just short of 7 
hours, the satellite will pass through the same geomagnetic 
longitude every 2 days. If the plasmasphere is moving close 
to the corotation speed, the same general region of the 
plasmasphere will be sampled at 2-day intervals. Therefore 
the observations in Figure 5 are assumed to represent four 
passes through two separate regions of the plasmasphere, 
one pair (Figures 5b and 5c) separated by 2 days and the 
other pair (Figures 5 a and 5 d) separated by 4 days. 

The recent history of magnetic activity, as measured by 
Kp, for the time period of these observations is shown in 
Figure 9. The Kp value had not been high for more than 10 
days prior to the first of these observations (reaching 4 on 
only a portion of November 23), and it had been steady, 
with daily sums near 24. On November 25, at the beginning 
of this sequence of passes, Kp increased from 2- to 5 in a 
period of 15 hours, then decreased and remained low 
throughout the remainder of these observations, with values 
less than 3+ and daily sums of 18 or less. The region in 
Figure 5 a was traversing the dayside on November 25 when 


the Kp began to increase and was in the nightside sector 
when the magnetic activity, and the inward convection, 
began to decrease again. Thus while there was loss of 
plasma in these tubes on the dayside, there was probably 
little compression of the plasmasphere in this region on 
the nightside, and the plasmapause remained at a relatively 
high L value of 5. The velocities in this pass show a 
variation with L similar to those in Figure 4c, with both 
ions having approximately the same speed at a given L value 
and the speed increasing with increasing L. In contrast, 
the region observed in Figure 56 was passing through the 
formative night sector when the magnetic activity began to 
increase, resulting in compression of the plasmasphere as 
evidenced by the location of the plasmapause at !=4.4. 
The velocities in this pass were higher than in Figure 5a for 
L< 3-8 but lower for !>3.8. While the flow in both cases 
appears to be from the southern hemisphere which is fully 
sunlit (as opposed to the northern hemisphere in which the 
shadow height was 200 km), the higher velocities for K3.8 
in Figure 5b may be indicative of an enhanced downward 
flow from the equatorial plasmasphere due to this com- 
pression [ e.g., Park, 1970]. 

The results in Figures 5c and 5 d represent subsequent 
observations of the regions of the plasmasphere first 
shown in Figures 5b and 5a , respectively. In the 4 days 
between the observations in Figure 5 a and those in Figure 
5d the velocities decreased throughout the region of 
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Fig. 9. The Kp history for the regions of the plasmasphere observed in Figure 5. The time lines at the top of the figure 
are computed as in Figure 8. 


observation and the plasmapause position moved outward 
from L- 5.1 to L- 5.5. Note that outside of L- 4.3 in Figure 
5d the He+ velocity is greater than the H + velocity. For 
the regions in Figures 56 and 5c the plasmapause position 
also moved outward with time from L-4.4 to 1=5.0 and 
between L= 4.2 and 1=5 in Figure 5c, H + and He + were 
observed to have significantly different velocities. The 
differences in the He + and H + velocities in the outer 
plasmasphere are indicative of the difference in the source 
strengths for the two ions. 

Substorm-Induced Effects (Refilling) 

The ion flux in the outer plasmasphere, where flux tubes 
are rarely full, is generally discussed in terms of the high- 
speed, often supersonic, flows associated with the early 
stages of refilling [e.g., Banks and Holzer , 1968; Nagai 
et al, 1984]. However, the plasma becomes subsonic in 
most of the flux tube within 2 days of the beginning of 
refilling [Banks et al., 1971; Khazanov et al, 1984 ,Horwitz 
et al , 1984] , and while the initial stages of refilling will be 
characterized by high-speed ion beams, the later stages (and 
indeed the dominant stage in terms of plasma supplied to 
the tubes) should be characterized by subsonic plasma with 
a nearly Maxwellian distribution f Horwitz et al., 1984]. 
Richards et al. [ 1 983 ] predicted that during this later stage 
of refilling the H + flux should be upward from both hemi- 
spheres into the equatorial reservoir, with He + flowing from 
the winter hemisphere, across the equatorial plane, down to 
low altitudes (<6000 km) in the summer hemisphere. An 
example of such refilling flows in the outer plasmasphere 
was shown in Figure 4b. In this case both H + and He + 
were observed flowing upward in the northern or fall 
hemisphere, consistent with the predictions of Richards 
et al [1983]. The magnitudes of the velocities are also in 
general agreement with Richards et al. 

Additional observations of high-speed flows, in what 
are believed to be refilling regions, have also been made. 


The results in Figure 6 are an example. These observations 
were made in the same general region of the plasmasphere 
with each observation separated in time by 2 days. This 
covers periods before an increase in magnetic activity 
(Figure 6 a), 1 day after the peak of this increased activity 
(Figure 6b) and 3 days after the peak (Figure 6c). Between 
the observations in Figures 6a and 6 b, this region of the 
plasmasphere passed through the nightside sector during 
two periods of enhanced magnetic activity - one on 
October 14 when Kp reached 8 and another on October 
15 when Kp reached 6- (see Figure 10). The high convec- 
tion electric fields associated with this activity would 
have compressed the plasmasphere in this region. This 
is evidenced by the change in the plasmapause location 
from L=4.5 in Figure 6 a to 1=3,5 in Figure 6b. Following 
the large storm on October 14, Kp began to decrease, 
finally reaching values less than 2 on October 16. Thus 
when observed on October 17 (Figure 6c), these flux tubes 
had been refilling for 2 full days. The evidence of this 
refilling is in the high H + speeds (up to 1 km s'l) in the 
region between the plasmapause location in Figure 6b 
and that in Figure 6c (/,= 5.0). 

The flux, in the refilling region in Figure 6c, at L- 4.5 
and an altitude of about 15,000 km was 2x1 0 7 ions 
cm' 2 s‘l for H+ ([H + ]=200 cm* 2 ). The values for He + 
were 5x1 0^ cm*- s - l and 80 cm’3, respectively. In Figure 
4b (at 1=3.5) the values are 10 cnr 2 and 2xl0 5 ions cm* 2 
s’ 1 for H + and 3 cm'3 lxlO 5 ions cm' 2 s~l for He+. The 
density values used in these comparisons are not estimated 
but have been taken from a complete analysis of the RIMS 
data (R. H. Comfort, personal communication, 1985). 

In the observations in Figure 6c, H + was observed flow- 
ing downward in the northern (autumn) hemisphere. If 
this flow is assumed to be interhemispheric, it differs in 
direction from the observations in Figure 4b in which H* 
and He + were observed moving upward from the northern 
hemisphere. The results in Figure 4b are more in line with 
the theoretical predictions of Richards et al. [1983]. The 
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Fig. 10. The Kp history for the regions of the plasmasphere observed in Figure 6. The time lines at the top of the 
figure arc computed as in Figure 8. 


region in Figure 4 b was in an earlier stage of refilling than 
the region in Figure 6c, as evidenced by the density of H + 
(1-10 cm'3 in Figure 4b and 100-1000 cm~3 in Figure 6c 
(R. H. Comfort, personal communication, 1985)). This 
may explain the difference in flow direction on the two 
days as well as the fact that in Figure 6c the maximum 
velocity was higher than in Figure 4b. Since the southern 
hemisphere was sunlit earlier, there should have been a 
tendency toward south-to-north flows. This would enhance 
the northward H + flow in Figure 6c while hindering the 
southward H* flow in Figure 4b. It should also be noted 
that in Figure 4b the He + velocity exceeds the H+ velocity, 
as predicted by Richards et al [1983]. It is also possible 
hat the observations shown in Figure 4b were made near 
:e flow reversal where steep gradients in the velocity are 
expected [Richards et al, 1983], Thus the observed 
velocity may be significantly lower than the maximum 
velocity along the field line. 

- He* Counterstreaming 

Richards et al, [1983] predicted counterstreaming of 
H + and He+ in the plasmasphere during refilling. The 
mechanisms responsible for this behavior in the flow of 
H + and He + are the winter anomaly and the winter helium 
bulge which produce H + density asymmetries between the 
winter and summer hemispheres of a factor of 2-3 and He + 
asymmetries of as much as an order of magnitude. The end 
“• tilt is that the partial pressure of all ions, as well as the 

al plasma pressure, is larger in the winter, and the winter 
litllow of ionization is able to push across the equatorial 
plane into the summer hemisphere before being offset by 
the summer hemisphere outflow. The effect is largest for 
He + which can lead to H + -He + counterstreaming in the 
summer hemisphere, as predicted by Richards et al 
[1983]. 

In the observations in Figure 6c, which were made 
during refilling conditions, the velocity values for He + , 
in the outer plasmasphere, are of the order of 100 m s'l 


(the quoted maximum expected error determined in section 
3). Nevertheless, they are consistent in magnitude and 
direction (from north to south), suggesting that H + and 
He + are indeed counterstreaming in this situation. In 
contrast to the results of Richards et al [1983], this 
counterstreaming occurred in the fall hemisphere. In this 
case the winter anomaly was probably not at its maximum 
(which occurs near solstice) so that the difference in the 
H+ source strength in the respective hemispheres may not 
have been as significant as it was in their model. Further- 
more, an effect which competes with the winter anomaly 
is that of the different solar zenith angles in the conjugate 
hemispheres. The stronger solar flux in the spring hemi- 
sphere results in an interhemispheric flow of H + which 
simulations show can still be significant at the local time 
(0930 hours) of these observations [cf. Bailey et al, 1978] . 
The extra winter enhancement for He + (the helium bulge) 
may result in He"*" being able to resume its winter to 
summer flow at earlier local times than H + , resulting in 
counterstreaming in the fall hemisphere with He + upward 
and H + downward as observed. 

Additional observations of counterstreaming are pre- 
sented in Figure 7 and show He + flowing upward in the 
winter hemisphere while H + is flowing downward (in one 
case, H+ is essentially stationary). These observations 
differ from those in Figure 6c in that they were made in 
the late afternoon to evening sector (1500-1900 LT). The 
two cases from November 25, 1981 (Figures la and lb) 
occurred in the postsunset period and show H + flowing 
downward (or south to north) and He + flowing upward 
(or north to south). These flows seem to be indicative of 
refilling but may also be influenced by the later sunset in 
the southern hemisphere, resulting in H + flow from that 
hemisphere. It is also plausible that the H + is actually 
flowing downward due to the collapse of the postsunset 
plasmasphere. 

In the case of January 25, 1982 (Figure 7c), prior to the 
observations there had been an extended period of more 
than 3 days during which Kp was approximately 4. In the 
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last 12 hours before the observations the Kp had dropped 
to 2. This suggests that these observations were made in 
a refilling region. Notice that this case is earlier in local time 
than the others and also is at lower altitudes in the winter 
hemisphere. As before, the counterstreaming is occurring 
in the winter hemisphere and thus is not in line with the 
theoretical results of Richards et al. [1983] for refilling 
situations. However, using a simple centered-tilted dipole 
magnetic field model shows that the difference in the 
solar zenith angles between the two hemispheres yields a 
factor of 2 larger solar flux in the southern ionosphere. 
This could overcome the larger H + production in the 
winter hemisphere due to the winter anomaly (also about 
a factor of 2) such that the resulting flow is from the 
summer to the winter hemisphere. It should also be noted 
that the observed velocities (maximums of 400 m s‘l for 
H + and 800 m s'l for He+) are an order of magnitude 
higher than those found by Richards et al [1983] whose 
simulations were for a flux tube which had been depleted 
to only one-half its steady state content. The higher 
observed velocities are more in line with the results of 
Khazanov et al [ 1984] for H + in a more severely depleted 
flux tube. Khazanov et al. did not include He + in their 
model, so no direct comparisons, with regard to counter- 
streaming, can be made. This difference in the stage of 
refilling could also be important in understanding why 
counterstreaming was observed in the winter hemisphere 
but not predicted in the Richards et ai. simulations. 

6. SUMMARY 

The observations of field-aligned ion velocities in the 
plasmasphere presented in this paper illustrate the dynamic 
nature of the plasmasphere under the influence of the 
diurnal cycle and varying levels of magnetic activity. 
Observations of large-scale variations in the shape and 
content of the plasmasphere induced by variations in mag- 
netic activity have been shown. The observations from late 
November 1981 (Figure 5) showed distinct differences in 
the ion velocities at different UT-LT sectors following 
enhanced magnetic activity. In addition, while velocities at 
both locations were observed to decrease as the plasma- 
sphere quieted and filled, these differences were retained 
over a period of 2-4 days. 

In most of these observations a summer-to-winter or 
spring-to-autumn pattern was evident in the flow of both 
ions, except in the cases of counterstreaming in which 
He + was moving from winter to summer. Also there was 
less variation with L in the velocities in the inner plasma- 
sphere (T< 3) than in the outer plasmasphere. This behavior 
is expected since the inner plasmasphere is less affected by 
magnetic activity and remains near equilibrium at most 
times. Typical speeds in the inner plasmasphere ranged 
from almost zero up to about 600 m s~l, while in the outer 
plasmasphere, speeds up to 1 km s'l were observed. The 
highest velocities (approaching 1 km s'l) were observed in 
regions of refilling. In one case (Figure 4 b) both H+ and 
He + were observed flowing out of the autumn hemisphere, 
consistent with the theoretical work of Richards et al 
[1983]. In the inner plasmasphere the observed velocities 
were somewhat higher than those predicted in simulations 
[e.g., Bailey et al, 1978]. 

Further observations of ion counterstreaming at high 


altitudes in the outer plasmasphere show H + flowing down- 
ward with He + flowing upward in the fall/winter hemi- 
sphere. While these observations appear to be from regions 
of refilling, they are not consistent with the theoretical 
work of Richards et al [ 1983] . The difference is attributed 
to differences in the solar zenith angles between the two 
hemispheres which tends to oppose the effects of the 
winter anomaly on H + . 

These results are significant both qualitatively and 
quantitatively for several reasons. First, they represent the 
first direct observations of the ion velocity at high altitudes 
in the plasmasphere. Second, these results illustrate the 
necessity of velocity measurements in conjunction with 
density, temperature, and composition measurements for 
understanding the dynamics of the plasmasphere. Velocity 
measurements in this context add an additional level of 
distinction between different states of the plasmasphere. 
Thus to fully understand the processes of plasma transport 
within the plasmasphere, the velocities must be measured. 
This points out the need for future instruments designed 
specifically to measure ion velocity in the plasmasphere. 

APPENDIX 

The flux of ions^/ at any location in space at time / 
is given by 

Jijj t ) = <«,•> / £ fjj.v, t ) d 3 v ( 1 ) 

where fJr,K t) is the distribution function. Writing this 
equation in terms of the speed and a solid angle gives 

J(r,t) = <n t > fdD.fv 2 dvfv (2) 

For the RIMS instrument the detector count rate can be 
written in terms of the distribution function, i.e., 

CR=<n l >fd£lfv 2 dvfvR(v) (3) 

Where R(v) is the velocity dependent instrument response 
function. For a cold ion population, ,R(v) can be taken to 
be constant, and the count rate can be approximated by 

CR = G <n t > f v 2 dvfv (4) 

where G is a geometrical factor which includes the response 
function and the effective solid angle of the instrument. 
For an integral detector such as RIMS, the integral over 
the magnitude of the velocity is done by the instrument 
in the zero retarding potential mode. In addition, since 
RIMS has an aperture width in the direction perpendicular 
to the spin plane of 40°-50° (FWHM), this portion of the 
total solid angle integration is also carried out by the instru- 
ment. The final integration is over 0, the angle in the spr 
plane of the spacecraft (see Figure 2), i.e., 

Jl %t) = f d<t>F v (<t>t] ( 5 ) 

where 

CR 

F v (<t>t) = -=jr = <■«,> / v* dv fv (6) 

is the quantity measured by RIMS. This yields 
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J(r,tJ = n t <^ tQt > 


(7) 


The Editor thanks C. Beghin and J. J. Sojka for their assistance 
in evaluating this paper. 


where £tot includes all naturally occurring velocity com- 
ponents as well as the satellite ram velocity. Returning to 
equation (3) and changing to components along the space- 
craft x and y axis yields two integrals, 


I x = / l£ y l cos <t> dtp 

( 8 ) 

ly = n i / !£ v l sin 0 dtp 

(9) 


From these two components the direction of ^tot can be 
found, i.e. , 


tan 0 Q 



<itot>y 




( 10 ) 


where 0 O is the angle from the spacecraft x axis to £ tot* 
Finally, the total velocity and the velocity along B can be 
found from the relation 


~B ~tot “~ExB ~ram 


(ID 


if Jiram an d iExB (which includes convection and corota- 
tion) are known. Separating ^ExB and ilrarn into ^rx ancl 
v ry gives 


V B X = Ji’tot' cos< ^o - v rx = b 1 cos b 


( 12 ) 


''^ = ^tot lsin ^o- V = b lsin b 


(13) 


where $3 is the angle of B with respect to the x axis. Solv- 
ing these two equations for^g gives 


b l = 


v - v tan 0 
ry rx v o 


(14) 


sin 0 g - tan 0 o cos 0 ^ 

be found from eit 
iity is defined by 

/ \F y \ cos (<H> 0 ) d<p 


iltot can then be found from either (10) or (11). An 
approximate density is defined by 


litot 1 

and the approximate flux along B is given by 


(15) 


$B = n £B 06 ) 

Because of the presence of a positive floating spacecraft 
potential, resulting in the exclusion from the instrument 
of some ions, this type of determination of both the 
density and flux leads to lower limits for both. 
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Update on polar wind thermal characterization scheme 

The application of thermal analysis techniques to the DE 1/RIMS polar wind dataset 
is an interesting challenge. While the RIMS design was a significant improvement to 
previous Retarding Potential Analyzer (RPA) instruments, its mass-resolving capabilities 
introduced additional complexities in the form of an energy-dependent solid angle aperture 
and a finite energy bandwidth beyond the RPA setting. In addition, in-flight anomalies 
produced significant difficulties in comparison to the analysis of conventional RPA 
datasets. Particularly relevant to the RIMS polar wind dataset is that the radial head RPA 
mechanism failed to operate properly early in the mission lifetime. 

Fortunately though, by adopting a numerical quadrature approach to the solution of 
.the response integral, using sums over the aperture area, solid angle, and energy limits, 
energy-dependent complications may be sufficiently simulated. And, because separate spin 
and RPA profiles are available for each of the three instrument heads, the radial head spin 
curve and the end head RPA curves are usable for temperature analysis even in the absence 
of the radial head RPA information. Since the effect of the radial RPA failure is as if the 
retarding voltage was constantly zero, the. total counts measurement remains accurate. 
Numerical algorithms of this type tend to be computationally intensive, hut have the 
advantage of being applicable to a variety of plasma distribution functions. 

In the case of the polar wind dataset, further considerations on the data analysis 
technique are necessary. Strictly speaking, spacecraft spin angle distributions observed in 
situ are not true pitch angle distributions. The measured spin angle is defined in the 
reference frame of the spacecraft, which is different from the plasma rc3t frame due to 
spacecraft motion (2 to 8 kilometers per second, depending on spacecraft altitude). This 
difference In reference frame would be of little consequence provided the thermal speed of 
the plasma is large compared to the spacecraft "ram speed. Fur a spacecraft median speed 
of 5 km/s, a faetor-of-five ratio of thermal speed to spacecraft speed implies 3.5 eV fur H + , 
13 eV' for He + and 52 eV for 0 + . Because polar wind energies are considerably below this 
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threshold, the RIMS spin angle data contain information about both the angular distribution 
of the plasma in its own frame of reference (the conventional pitch angle distribution; and 
the drift motion of the plasma relative to the spacecraft. These two aspects of instrument 
response tend to dominate independently - one at times obscuring the other - depending on 
the plasma's thermal speed and its relation tc the spacecraft and plasma drift speed. 

Roughly speaking, it is for subsonic (thermal speed greater than drift speed) ion 
Mach numbers that the angular characteristics of the distribution arc best resolved and for 
supersonic Mach numbers that drift motion is most easily resolvable. Note that in the case 
of H+, flow will be subsonic for thermal components as low as 1 eV when drift velocities 
are less than 15 knVs. The bottom line, for KIMS polar wind data, is that the relative drift 
motion will dominate the spin angle response, leading to a "rammed" flux peak in tne 
direction of relative wind as seen in the spacecraft frame. The angular (and hence, 
temperature) information will be apparent only in the detailed shape of the rammed peak 
and may, at times, become indeterminate. The end-head RPA profiles then, become 
increasingly important to the accurate determination of pur polar wind thermal properties. 

With all these, cautions in mind, a numerical simulation model of the RIMS 
response to a flowing plasma distribution was written [Giles etal , 1 99.1]. The simulation 
technique improves on previous work by realistically modeling the interaction of a moving 
spacecraft with a drifting plasma and by allowing insertion of a variety of distribution 
functions (for example, bi-Maxwcliian, generalized Lorentz (Kappa), end Spitzer-Hdrm 
[1953]), the three ‘Component drift velocity, and a variable spacecraft potential. For the 
polar wind study, the distribution function is assumed to be a flowing bi-Max weiliar. which 
Is accelerated through the potential difference between the spacecraft and the plasma. The 
potential drop 13 assumed to occur in an infinitesimally Ulin sheath about the spacecraft and 
the detector, rendering the detail of sheath structure unimportant. The implementation 
allows simultaneous fitting of the radial sector spin curve data and the end see, or RPA 
(energy) curves. 
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Figure 1 Rhows simulated RTMS spin-angle and RPA responses lo a bi-M&xwellian 
plasma distribution (solid lines) and the RIMS polar wind data (crosses) for 1982/106 
22: 17 UT. Note that the instrument response, for a given detector look direction, is to ions 
with velocities anti-parallel to the detector look direction Plasma distribution parameters 


for this fit arc: 


density 

224 cm 3 

T|| 

0.36 eV 

T± 

0.28 eV 

field-aligned velocity 

4.14 km/s 

anti-sunward convection velocity 

2km/s 


Note that these parameters correspond to a cold subsonic outward flow (ion Mach 
number -0.49). Early predictions indicated that the H + polar wind thermal distributions 
should be anisotropic, with Th>Tj. for supersonic outflows and Th<Tj. for subsonic flows. 
The present example however, fits in with work by Wu and Taieb [1993] showing that an 
upward heat flux may exist even for subsonic flow at low altitude in the high-latitude 
ionosphere. 

Processing of the. RTMS polar wind dataset using the RIMS simulation code 
continues as of April, 1995. Several studies have been completed which relieve difficulties 
associated with fitting multivariate functions. Frequently it is a matter of art to decide 
which variables in a model should be fit from the dataset, and which should be held 
constant at fixed values, for example values predicted by theory or derived in previous 
analysis, In the case of the polar wind, data set, variables potentially free for fitting to the 
data are the spacecraft potential, plasma density, the two-componcnt plasma temperature, 
and the three-component plasma bulk flow velocity. In the present implementation, the ion 
field-aligned velocity component and flux measurements (and so, the density, value) are 
held fixed. Their values were obtained [ Chandler et al, 1991] by integrating the observed 
flux over the spin cycle of the satellite for a zero RPA voltage. Alsu, the spacecraft 
potential value is held fixed during the first stage of the fitting process, an initial value 
being derived from an empirical relationship between the plasma density and spacecraft 
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potential [ Comfort et al, 1988], Values fit from the data set then are the two-component 
temperature and the remaining two components of the bulk flow velocity. The 
perpendicular temperature is fit by chi- square minimization to the end-head RPA profiles. 
The parallel temperature is similarly fit to the radial head spin angle profile. The orientation 
of the spacecraft for the polar wind data set is such that the expected anti-sunsvard 
convection velocity is almost completely described by the y component of tire bulk flew 
velocity vector. This is fit first by minimization to the taJiul head spin angle pioflle and 
then refined through minimization to the end-head profiles. Last, the spacecraft potential 
value is fit in the same order, to the same three profiles. The fit in figure 1 was obtained by 
this method. 
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